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Abstract 

Electrostatic plasma confinement and magnetic electro- 
static plasma confinement (MEPC) have been studied 
for four decades. The multiple potential well hypothesis, 
postulated to explain high neutron yields from Hirsch’s 
colliding beam experiment, has been supported by sev- 
eral pieces of evidence, but results were inconclusive. 
Magnetic shielding of the grid was developed to reduce 
the required beam current and to prevent grid overheat- 
ing. Electrostatic plugging of magnetic cusps evolved to 
a similar configuration. Due to low budgets, early 
MEPC experiments used spindle cusps, which are poor 
for plasma confinement. Later experiments used multi- 
pole cusps or a linear set of ring cusps, which have larger 
volumes of field-free plasma. To keep the self-shielding 
voltage drop A 0 < 100 kV, the electron density in the 
anode gap should be less than about io^nD. The central 
plasma density can be an order of magnitude higher. The 
ATOLL toroidal quadrupole had anomalous electron 
energy transport, but the JupiteraM linear set of ring 
cusps achieved a transport rate about a factor of two 
above the classical rate. With near-classical transport, a 
power gain ratio Q = xo is predicted for a reactor with r p 
= 3 m, B = 6 T, and applied voltage 0 a = qookV. Besides 
producing electricity and synthetic fuels, MEPC reac- 
tors could be used for heavy ion beams sources and neu- 
tron generators. The main issues of concern for MEPC 
reactor development are electron transport, plasma pu- 
rity and electrode alignment and voltage holding. 


List of Parameters 

(SI Units are used, unless otherwise noted.) 
a = half-width of the anode gap 

A = area 

A / oss = effective area for plasma loss out of magnetic 
cusps 

A = magnitude of the magnetic vector potential 
B = magnetic field (magnetic induction) 

Ba = magnetic field in the anode gap 

Bi, = [zp 0 nkfT^ + Ti)] lh = boundary magnetic field 

2U = peak magnetic field at the coil 

B ^ = point cusp magnetic field 

c = speed of light 

c s = ion acoustic speed 

d - grid wire diameter 

D a = electron diffusion coefficient due to collisions 
with neutral atoms 

D e i - electron diffusion coefficient due to collisions 
with ions 

D max = maximum diffusion coefficient 

D 0 = total collisional diffusion coefficient 

e ^ = electronic charge 

fi = fraction of incident neutrals which are ionized 
(instead of charge exchanged) 

fit, = fraction of fusion product alpha energy W”, 
transferred to electrons 

fai = fraction of fusion product alpha energy W”, 
transferred to ions 

F = geometrical factor having dimensions of 
reciprocal length 

g = spacing between grid wire centers 
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h = width of the anode along the magnetic field 

direction 

I = electron current from electrodes into the 

plasma 

Idiff = electron diffusion current 

It, = electron injection current 

= wave number = znfk 
= Boltzmann constant 

k „ = subscript denoting either electrons or ions 

L = length between ends of a linear cusp device 
Lx = distance between adjacent ring cusps 
m. , = electron mass 

tru = number of pairs of multipole conductors 

= N /2 

M = ion mass 

ti~> = plasma electron density 

n^a = uncompensated electron density at centre of 
anode gap 

tl>b = Brillouin density 

= particle density of species k 
ru c = central core electron density 

JVe, = total number of electrons 

tL>o = central plasma density 

71 -jo = initial plasma density (pulsed electrostatic 

confinement) 

(l,j = ion density in the anode gap 

= trapped cold electron density in anode region 
N = number of multipole cusp gaps 

p = ( T ±1 T=) l/2 

Pe = electron heating by the edge electric field that 
affects the central plasma it = ion heating by 
plasma waves 

Prod = radiative power loss 

P-w = electron power loss to plasma waves 

P 0 = canonical angular momentum 

Q = (fusion power)/(input power) 

q = C0 2 p e /0l) 2 ce 

r = radius 

r a = point cusp anode radius 

n = plasma boundary radius 


r„ = virtual electrode radius (electrostatic con- 
finement and Penning trap) 

r 0 = radius of plasma boundary surface in point 
cusp 

r p = approximate plasma radius at magnetic field 
boundary 

R = accelerating grid radius (electrostatic 
confinement) 

R = magnetic mirror ratio 

R = ring cusp anode radius 

S = plasma surface area 

Si = ionization source 

Sm, = volume-averaged neutral input rate 

T e ^ = electron temperature (keV) 

TL = electron temperature 

T eo = central electron temperature 

u = ion speed 

Uia = ion speed in anode gap 

= electron velocity 

O’Va, = electron density-gradient drift velocity 

Y = plasma volume 

r w = characteristic half-width of the electron 

distribution in anode region 

W a = average energy acquired along the slope of the 
potential well by ions 

Wt, = electron energy 

Wj = ion energy 

Wtoss = average electron energy loss per incident 
neutral atom due to ionization 

W_l = ion energy perpendicular to the radial 
direction 

x -r/R 

y = OJTt, 

= zg,(DkfI \ 

Y = inR/N 

Z\ = charge number 

Zk = x for electrons and 1/2 for ions 

a = a parameter depending on T/ZV, 

Y = fraction of current intercepted by the grid 
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5 

= the effective half-width of the cusp gap for 
plasma loss 

%i 

= ion confinement time 


To 

= electron collision time 

5 

= a function, of Ph P2 

x trap 

= time for central plasma free electrons to be 

A 0 

= self-shielding voltage drop in anode gap 


trapped by edge magnetic field 

A r 

= thickness of plasma edge region 

<t >A 

= applied voltage 

e 

= (trapped electron pump-out rate)/ 


= potential barrier for electrons or ions 


(diffusion rate) 

<ZL 

= electron potential barrier 

8 

= effective open fraction of wire mesh grid 

0 i 

= ion potential barrier 

8e 

= geometrical open fraction of wire mesh grid 

0 P 

= plasma potential = 0 ,- + A 0 

£0 

= permittivity of free space 



1 

= global Brillouin ratio 

(0„ 

= cyclotron frequency 


C Oce, 

= electron cyclotron frequency 

In A 

= Coulomb logarithm 

COn 

= ion cyclotron frequency 

l 

= wavelength 


P» 

= permeability of free space 

0 )pe^ 

= electron plasma frequency 


= electron-ion momentum-transfer collision 




frequency 



P* 

= electron Larmor radius in anode gap 



P* 

= electron Larmor radius at plasma boundary 



P^ 

= electron Larmor radius 



P» 

= electron Larmor radius in point cusp 



P* 

= ion Larmor radius 



P* 

= inner electron turning radius 



P* 

= outer electron turning radius 



<O r V> 

= DT fusion reaction rate parameter 



T± 

= electron loss time by diffusion across the 




magnetic field 



T 

= pulse length 



Xk 

= characteristic collision time for electrons or 




ions 



w 

= time for heat loss by conduction 



W 

= characteristic time for electron diffusion 




across the magnetic field 




= electron lifetime or confinement time 



TJS± 

= electron perpendicular energy confinement 




time 



Teh 

= electron parallel energy confinement time 



T e i 

= electron-ion collision time 



T eq 

= electron-ion equipartition time 
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i. Introduction 

Magnetic electrostatic plasma confinement (MEPC) 
uses electrostatic fields to enhance plasma confinement 
by open magnetic fields. The purpose of this review is to 
describe: 

(0 the evolution of MEPC from electrostatic 
plasma confinement and from open magnetic 
confinement systems 

(ii) MEPC theory 

(iii) experimental observations 

(iv) fusion reactor concepts 

(v) salient issues 

MEPC devices have also been called ‘electromagnetic 
traps’, ‘magneto-electrostatic traps’, ‘magnetic inertial- 
electrostatic confinement’, ‘electrostatically plugged 
cusps’, ‘kinetic-electric-magnetic plasma confinement 
(KEMP)’, ‘electric-magnetic confinement’, ‘Polywell’ and 
‘modified Penning trap.’ Discussions of electrostatic 
field applications in tandem mirrors, bumpy tori and 
tokamak helicity injection are beyond the scope of this 
review. 

1.1 Electrostatic Plasma 
Confinement 



Deflections of injected particles towards nearby grid 
wires decrease the effective open area of the grid by an 
amount 7 


£ 

£ e 



(x) 


1.1.1 Virtual Electrode Formation 

Lavrent’ev 1 first proposed electrostatic plasma confine- 
ment on 22 June 1950, and magnetic field enhancement 
of electrostatic plasma confinement in March 1951. In 
the early 19605 Farnsworth 2 proposed electrostatic 
plasma confinement for the production of fusion reac- 
tions, based upon his experience with spherical multi- 
pactor vacuum tubes. Lavrent’evx developed a theory of 
charged particle flow and focusing in plane and spherical 
geometries. Budker4 calculated the fusion reaction rate 
that could be attained by trapping ions in the negative 
electrostatic charge of an electron beam. Early electro- 
static plasma confinement experiments^ consisted of 
concentric spherical wire mesh electrodes, which accel- 
erated charged particles inwards to produce a virtual 
electrode inside the central sphere (Figure 1). Charged 
particles of the opposite sign would be attracted into 
the virtual electrode. For example, if electrons were in- 
jected into the sphere, they would produce a virtual 
cathode 6 , which would attract positive ions. In order to 
minimize grid current and heating, it is desirable to have 
the grid mesh nearly transparent to the charged particle 
flows. 


where e is the effective open fraction of the grid, e g = (g - 
d) 2 /g 2 is the geometrical open fraction of the grid, g is the 
spacing between grid wire centers (square mesh) and d is 
the grid wire diameter. For example, ifg = 4 mm and d = 
0.2 mm, then z g = 0.90 and e = 0,865. The circulating 
current I c is related to the current I intercepted by the 
grid wires by the equation: 


For the example case with 8 = 0.865, the theoretical 
value of IJI = 6.9. Thus, for a highly open grid the circu- 
lating current can be about an order of magnitude 
greater than the grid current. 

Lavrent’ev and co-workers 8 studied ion generation by 
oscillating electrons and ion focusing with spherical 
grids. Probe measurements inside a 5 cm radius spherical 
cathode indicated a strong focusing of injected ions 
down to radii < 1 cm. The focused ion current density 
was 12-25 mA cm' 2 at 4 kV in 0.13 Pa of argon . 9 f° A 
passing ion density of io J 7 m' 3 ; was deduced from cur- 
rent probe data. 12 
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Hirsch'? injected electrons into a spherical wire mesh 
anode, and used an electron beam probe to measure the 
potential well depth produced by the virtual cathode. 

He found that deep, negative potential wells were easily 
produced and maintained. He also observed some oscil- 
lations of the discharge at frequencies ~ioo MHz. 

For the case of electron injection, to produce a virtual 
cathode, there are two great difficulties: the very high 
currents required 1 ^ and grid wire overheating. Better 
confinement could be obtained for the case of ion injec- 
tion into the accelerating grid, producing a virtual an- 
ode. 

Rockney '5 simulated the formation of a virtual electrode 
in cylindrical geometry by solving Poisson’s equation for 
a case with only ions, represented by 2000 charged rods, 
and showed that the virtual electrode is stable for long 
periods of time. In a simulation with both electrons and 
ions, however, Barnes 16 found that electrons trapped in 
the virtual anode tend to wash out the anode potential, 
so that the virtual anode gradually disappears. Only in 
cases with very high ion currents did the virtual anode 
persist. He concluded that the electrostatic-inertial 
confinement device, as originally conceived, was unlikely 
to confine a fusion plasma adequately Porter and Kie- 
vans 1 ? investigated the stability of electrons flowing 
among ions within a virtual anode, and concluded that 
the electrons are marginally, stable to the two-stream 
instability 

1 .1 .2 Multiple Well Hypothesis 

Hirsch 18 observed high DT neutron yields (up to 2 x 10 10 
s' 1 ) from a steady-state device using ion injection from 
six ion guns. These neutron yields exceeded the predic- 
tions of simple theories. One hypothesis to account for 
the observed neutron yields is the formation of multiple, 
concentric spherical potential wells of alternate sign 
inside the central sphere, with much higher circulating 
currents in the inside layers, as illustrated in Figure 2. 
According to such a model, ions trapped in interior lay- 
ers could have high energies, circulating currents and 
densities, and could produce significant fusion energy 
yields. Hu and Kievans 1 ? modeled multiple wells for 
cases with zero angular momentum (perfectly radial mo- 
tion). They calculated the attainable current amplifica- 
tion ratio (the ion current inside virtual anode/ion beam 
current outside virtual anode), and found that values as 
high as 150 were possible if large trapped ion popula- 
tions could be confined in the outer region. 




The following observations appear to be consistent with 
the multiple potential well hypothesis: 

(i) High neutron yield; 

(ii) Peaks in the curve of neutron yield versus radius 
measured by Hirsch 18 as shown in Figure 3; 
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Radial Position 




(deuterium)IS 


(iii) Peaks in the curve of electron bremsstrahlung 
yield versus radius measured by Hirsch' 8 as 
shown in Figure 4; 



-12.5 -6 0.0 6 12.5 19 

Position With Respect to Cathode Vertical Axis, mm 


(iv) A large time-delayed puff of gas after beam in- 
jection cut-off measured by Hirsch 18 ; 

(v) Theoretical studies by Cherrington et al 2 ° for 
cases with injected beams having narrow 
spreads of energy and angular momentum show 
a slight double well, as illustrated in Figure 5; 



Figure 5 — Theoretical variation of dimensionless electrostatic potential Y versus dimen- 

parency = 90%, and very small spreads of energy and angular momentum. A double well 
is apparent. 

(vi) Electron beam probe studies by Swanson et, 
al 2 '’ 22 of a spherical electrostatic confinement 
device at low pressures; 

(vii) Peaks in the curve of electron density versus 
radius (Figure 6) measured by laser heterodyne 
studies in a cylindrical electrostatic confine- 
ment device by Meeker and co-workers. 2 }-^ This 
figure is for the case of electron injection; simi- 
lar, less pronounced peaks were observed during 
ion injection. 
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erodyne interferometer during electron injection into a cylindrical grid with radius 40 mm, 

Other studies, however, have yielded results less indica- 
tive of multiple potential wells. Imel 2 5 and Black and 
Kievans 26 developed theoretical models of the electro- 
static potential profile and electron and ion density pro- 
files, using distribution functions tailored to various ex- 
perimental conditions. They analyzed three cases 2 7 

(i) Ion beam injection at very low pressures (the 
Hirsch experiment). They concluded that mul- 
tiple deep potential wells would probably not 
form under the conditions of the Hirsch ex- 
periment. 

(ii) Ion injection at higher pressures, where charge 
exchange produces a broad spread of ion ener- 
gies (several experiments). They predicted a 
central plasma density close to the value meas- 
ured in the Penn State experiment (2 x io 1 * m3), 
and found that substantial potential wells could 
be formed for p < 0.13 Pa (1 mTorr). 

(iii) Electron injection at low pressures (the Swan- 
son etal and Black and Robinson experiments). 
They predicted deep potential wells for high 
perveance cases, as observed by Swanson etal, 
and shallow potential wells for low perveance 
cases, as observed by Black and Robinson. 28 

They did not report any cases with finite angular mo- 
mentum and multiple potential wells. 

Gardner and colleagues 2 ? continued the ion beam ex- 
periment of Hirsch and used a microwave cavity to 
measure the number of free electrons in the plasma, but 


they were unable to reproduce the neutron yields ob- 
served by Hirsch. Baxter and Stewart3° predicted neu- 
tron yields close to those observed in the Hirsch ex- 
periment by analyzing ionization and charge exchange in 
the ion beam, without invoking the multiple potential 
well hypothesis. 

Peterson and Oleson? 1 studied hemispherical focusing of 
a low-energy ion beam in a double-beam plasma device. 
When the ion beam energy exceeded T c the beam ions 
were slowed down by a potential hill (virtual anode) at 
the focus. They observed formation of a dense group of 
cold electrons (0.8 eV), which appeared to ride along 
with the ion beam, in addition to the background 
plasma electrons (4.5 eV) and primary electrons from 
filaments. Such cold electrons could tend to wash out 
potential wells. 

Nadler et all 1 built a collimated solid state proton detec- 
tor to observe protons resulting from DD fusion reac- 
tions in a spherical electrostatic confinement device. 

The energy shifts produced by aluminum and lead foils 
verified that they were detecting 3 MeV protons. At 12 
mA cathode current, 30 kV voltage, in 0.5 P a of deute- 
rium fusion reactions between the focused ion beam and 
the background plasma were dominant. In future ex- 
periments at higher currents, the higher proton count 
rates should permit improved angular resolution, which 
could facilitate the detection of multiple potential wells. 

1.1.3 Regimes of Operation 

Miley et a fa distinguished three regimes of ion injection 
into an inertial electrostatic confinement (IEC) device: 

(1) A discharge mode. In this mode a self-sus- 
taining plasma discharge (like a glow discharge) 
generates ions throughout the chamber volume. 
Free electrons are produced by ionization and 
by secondary emission from grids. This mode 
requires comparatively high gas pressures (>i 
Pa), which result in substantial charge exchange 
rates. 

(2) An emitter mode, based on a patent by 
Hirsch and Meeks 34 In this mode primary elec- 
trons from thermionic cathodes orbit in and 
out of an intermediate grid, ionizing back- 
ground gas, at pressures typically 0.1-1 Pa. Ions 
produced by ionization are then accelerated 
into the central cathode. This mode was used in 
experiments by Meeker and co-workers 2 4>35 to 
measure density profiles. Using this mode with 
deuterium gas, Miley et a fa observed steady- 
state neutron yields of io5 neutrons per second. 

(3) A beam injection mode. High-energy ion 
beams are focused into the central region, to 
produce a high-density virtual anode. This 


Plasma Physics and Controlled Fusion 36 (1994), pp. 1539-1593, doi:io. 1088/0741-3335/36/10/001 


7 


mode, which can operate at the lowest pres- 
sures, was used by Hirsch to obtain high neu- 
tron yields. This mode is most suitable for fu- 
sion reactors, because the low pressure mini- 
mizes ion energy loss by charge exchange. 


1 .1 .4 Reactor Prospects 

For steady state . > operation with radiative cooling of a 
tungsten grid intercepting i% of the ion current and 
focusing the ions to a virtual anode at radius r 0 = o.oi R, 
Lavrent’evs 6 found that a grid radius J?>iom would be 
needed in order to achieve a net energy output. For the 
case of a pulsed reactor the required grid radius is given 
by 


R > O.OOlyv^l 


ft 


/4 d® 


where y is the fraction of current intercepted by the 
grid, Tp is the pulse length, d is the grid wire diameter 
and 0.4 is the applied voltage (kV). For a case_ with y = 
o.oi, r„/R = o.oi, ®a = ioo kV, x p = o.i and d = o.ooi m, 
this equation gives R > 2.5 m. 

It is also possible to focus the injected ions temporally 
by increasing the voltage during injection. 37 Ions starting 
out later are accelerated more and catch up with the 
early ions, so that the ion density at the focus is much 
higher than the steady-state value for continuous injec- 
tion. For example, starting with background plasma at 
radii between S = im and 2 m with density n ^ 0 = 10 18 m"3, 
the focused density would be ru = 6 x io 2 5 m i. The con- 
verging ion beams could also be used to compress and 
ignite inertial confinement fusion target pellets. 

Although purely electrostatic plasma confinement might 
be attractive for a fusion reactor with a large radius or 
short pulse length, overheating of the grid wires is a se- 
rious concern. One way to reduce the heat load from 
charged particle bombardment of the grid wires is to 
shield them magnetically by passing high currents 
through them, as illustrated in Figure 7. If electrons are 
injected between the grid wires, very few will contact 
the grid wires directly Some electrons will gradually be- 
come trapped in the magnetic field and diffuse to the 
grid wires, but this process is orders of magnitude 
slower than direct bombardment. This magnetic shield- 
ing of grid wires, here called MEPC, has also evolved 
from open magnetic confinement systems. 



1 .2 Open Magnetic Confinement 
Systems 

Minimum/! magnetic confinement systems offer good 
MHD stability, but suffer from rapid plasma loss along 
magnetic field lines.3 8 Confinement times in simple 
magnetic mirror cells are typically of the order of ion- 
ion scattering times, and confinement times in magnetic 
cusps are tens of transit times, which typically makes 

q _ (fusion power) ^ ^ 

(input power) ^ 

for ordinary mirrors and cusps. In order to attain Q > 10 
without requiring very large size or extremely high mag- 
netic induction, it is desirable to find a means for plug- 
ging end losses from magnetic mirrors and cusps. Many 
means of Q-enhancement for open magnetic confine- 
ment systems have been studied, including tandem mir- 
rors, multiple mirrors, field reversed mirrors, Tormac, 
Surmac, plasma rotation, gas dynamic traps, radio fre- 
quency plugging and electrostatic plugging.394° Axi- 
symmetric cusps have been suggested as electrostatic 
end plugs for tandem mirrors .+m 2 43 

Several cusped magnetic fields have been considered for 
plasma confinement,^ as illustrated in Figure 8. The 
toroidal multipole shown in Figure 8 has N =6 cusps 
(hexapole). In addition to these configurations, arrays of 
point cusps on the surfaces of polyhedrons have been 
studied by Keller and Jones45 and by Sadowsk i . 4 6 47 48 .49 5° 
All of these cusp geometries have been considered for 
electrostatically plugged devices. 


Plasma Physics and Controlled Fusion 36 (1994), pp. 1539-1593, 


doi:io. 1088/0741-3335/36/10/001 



For plasma confinement by magnetic cusp fields, the 
loss area through a ring cusp is approximately nzRiib), 
where R is the ring cusp radius and 5 the effective half- 
width of the cusp gap through which plasma is lost. In a 
spindle cusp, the loss area through the point cusps is 
comparable, so the total loss area A/ ass = HnRb. The ion 
flux into the loss regions is approximately 0.25ml, where 
u is the ion speed, so the confinement time in an un- 
plugged magnetic spindle cusp device is approximately 

nV V 

T ~ = ( 5 ) 

0.25 nuA loss 2nuR8 

where V is the plasma volume. This time is typically tens 
of bounce times back and forth inside the central field- 
free region. Spindle cusps have very small plasma vol- 
umes. If 6 - pi then unfeasible large R and B would be 
required for a reactor. If 6 — (ptp,)' 12 (a hybrid gyro- 
radius), as indicated in some cusp experiments and in a 
computer simulation, ? I >? 2 ’?3 then a reactor using a ‘picket 
fence’ cusp (such as the spherical cusp and linear set of 
ring cusps in Figure 8) is feasible, but high magnetic 
fields and large radii are needed. 5 + 

For a linear set of ring cusps (Figure 8) with distance L, 
between the ring cusps, the plasma volume between 


cusps is V = nrp 2 L, where r p is the effective plasma radius, 
and the loss area Ai oss = 4Jtr8. For this case 


fL, 

T = ^ 


This equation, which estimates the plasma loss time by 
free-flow-out ring cusps with an effective loss width 28 
in each ring cusp, will be used later to estimate the effec- 
tiveness of electrostatic plugging. 

Lavrent’ev and colleagues?”*? proposed electrostatic 
plugging of cusps to reduce the loss rate of plasma flow- 
ing along magnetic field lines out of the cusp gaps, and 
conducted early experimental studies of this technique. 
If the plasma potential is highly negative relative to the 
walls, the ions will be confined electrostatically Only 
electrons will be energetically able to pass through the 
narrow cusp gaps (plus a few ions in the high-energy tail 
of the Maxwellian), so the width of the untrapped elec- 
tron stream flowing through the cusp gap may be re- 
duced to 8 - pt,. Cathodes outside the cusp gap reflect 
escaping electrons back into the plasma. Such electro- 
statically plugged cusps are equivalent to magnetic 
shielding of the grid wires of an electrostatic-inertial 
confinement device. In electrostatically plugged cusps, 
the ions no longer flow freely out of orifices, so the free- 
flow confinement time estimates do not apply Instead, 
the ions are confined electrostatically by the potential 
well of electron space charge, and electron confinement 
is limited by diffusion across the magnetic field. 

Moir discussed the use of electrodes to contain warm 
plasma for stabilizing open magnetic confinement sys- 
tems, such as a toroidal quadrupole cusp.? 6 The augmen- 
tation of toroidal magnetic confinement by strong elec- 
trostatic fields was studied by Daugherty etah i and by 
Stix. 5 8 -? 9 . 6 ° Jones 6 '- 62 suggested the use of electron or ion 
injection to create radial electrostatic potential varia- 
tions, instead of using plugging electrodes. 
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2. Theory 


2.1 Electrostatic Potential 

Magnetic electrostatic plasma confinement, which may 
be accomplished with any magnetic cusp configuration, 
is illustrated in Figure 9 for the case of a toroidal quad- 
rupole cusp. High-voltage electrode rings are placed in 
each of the cusp gaps, as shown. The anodes are biased 
positive, and the cathodes negative. Operation of the 
confinement system depends mainly upon the voltage 
<t>A applied between the cathode and anode, and little 
upon the relative location of the ‘ground’ potential. 
From Figure 9 it is evident that 

(t> A =(t> e + <t>i + (7) 


This potential shape is the same as that of a simple tan- 
dem mirror. 



Figure 9 — An electrostatically plugged toroidal quadrupole magnetic cusp and the resul- 
tant distribution of electrostatic potential ®. The broken line represents the potential with 


In a vacuum, the interior of the device will be near the 
anode potential (broken curve of Figure 9). Plasma may 
be produced by electron beam injection into low- 
pressure gas, by rf (radio-frequency) heating, by plasma 


gun injection, by laser-pellet heating, by neutral beam 
injection or by other means. Electrons exiting the cusps 
are reflected by the negative voltage of the cathode, but 
the exiting ions are not initially confined. Since elec- 
trons are lost less rapidly than ions, the plasma develops 
a slight negative charge, and the plasma potential be- 
comes negative relative to the anodes (the full curve of 
Figure 9). If the anode gaps are narrow enough that their 
potential is not entirely shielded out by the plasma in 
the anode regions, then the potential there will be near 
the anode potential and higher than the plasma poten- 
tial, forming a potential hill 0,- for the remaining ions. 
Only a slight fractional charge imbalance is required to 
set up a potential hill many kV high. Then only ions 
with kinetic energies greater than eJX>, can escape along 
the magnetic field fines out of the cusps. The two- 
dimensional potential distribution (Figure 10) is saddle- 
shaped in the anode regions, with the saddle point lying 
an amount A 0 below the anode potential. 


0 

-1 

-2kV 

-3 

-4 

-5 


Figure 1 0 - Hypothetical equipotential surfaces in the anode region, showing saddle, 



The potential hill for ions approaching the walls (at an- 
ode potential) is even higher, so most ions will go out 
through the cusp gaps as soon as they acquire enough 
energy to overcome the barrier 0,; no ions acquire 
enough energy collisionally to surmount the barrier (0, 
A 0 ) and reach the walls. The ions are electrostatically 
confined in a negative electrostatic potential well pro- 
duced by a slight charge imbalance. The ion Larmor ra- 
dius is less significant than usual, because the ions are 
electrostatically reflected in the boundary layer. The 
maintenance of a potential difference along the mag- 
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netic field lines is facilitated by the boundary condition 
that the walls are very close to the plasma only in the 
anode gaps. Some energetic ions are replaced by cold 
ions as a result of charge exchange with neutral gas. 

The central plasma region is uniform in density and 
temperature, free of electric and magnetic fields, and 
surrounded by a thin edge layer. The outside boundary 
of the edge layer is the magnetic surface along which 
electron orbits graze the anodes, and the inside bound- 
ary is the uniform central plasma region. The boundary 
layer thickness at any location with magnetic induction 
B is found by magnetic flux conservation to be 


A r~a\ 


( 8 ) 


where a is the anode gap half-width, R is the ring cusp 
anode radius and r the radius of the boundary layer. If B a 
= 5T, B = iT, R = 4 m, r = 3 m and a = 1.5 mm then A r = 10 
mm. This thickness is much less than the plasma size, so 
the plasma has a ‘sharp boundary’ 


2a 



lines. Dots represent ions. Electron trajectories are as follows: 1 = reflected off plasma 


boundary, 2 = passing through anodes and reflected by cathode, 3 = magnetically trapped 
in boundary layer, 4 = electrostatically trapped by anode potential well. 


2.2 Electron Orbits 

There are several types of electron orbits, illustrated in 
Figure ii: 

(1) Electrons that stay in the central plasma and are 
reflected geometrically from the convex edge 
layer. 

(2) Electrons from the central plasma that can pass 
through the anode gap regions and be reflected 
by the cathodes. 

(3) Electrons that are magnetically trapped in the 
edge layer. (Many of these are mirrorreflected 
away from the cusp gaps.) 

(4) Cooler electrons that are electrostatically 
trapped by the anode potential. Most of these 
electrons are produced by ionization of incident 
neutral gas atoms along the slope of the poten- 
tial well. These trapped electrons are deleteri- 
ous, because they increase the self-shielding 
voltage drop A<t>. 

Electrons may be produced by thermionic emission 
from an electron gun, by secondary emission from the 
plugging cathodes and by ionization of neutral gas. 

Much of the injected electron beam is reflected by the 
plasma and lost to the cathodes. This portion of the in- 
jected beam does not contribute to plasma density 
build-up. Yushmanov^ found that only a small fraction 
of the injected beam current is effectively trapped. The 
time required for plasma density build-up by electron 
beam injection (many particle confinement times) ex- 
ceeds the duration of the flat-top magnetic field in some 
pulsed experiments. 

Many of the electrons produced by ionization will be 
magnetically trapped in the edge region. The untrapped 
electrons will circulate back and forth through the 
plasma, occasionally passing out through the anodes and 
being reflected by the potential barrier <ty.. Coulomb 
collisions and electric field fluctuations cause diffusion 
in velocity space. Trapped electrons can become heated 
and detrapped, and free electrons can be deflected and 
trapped either magnetically or electrostatically. 6 * 

There are two main electron loss processes: trapping by 
the magnetic field with subsequent diffusion to the an- 
odes, and collisional diffusion in velocity space over the 
potential barrier with loss to the cathodes. (Loss to 
the cathodes converts electron kinetic energy into elec- 
tricity, returning energy to the power supply) Recombi- 
nation with ions is negligible. 
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2.3 Loss Rates over 
Potential Barriers 


The loss rates along the magnetic field by velocity-space 
diffusion over the potential barriers have been calculated 
by Cohen etal, 6 $ based on the method of Pastukhov 66 
For a plasma with singly charged ions, the resulting par- 
ticle^ loss rates along the magnetic field are 

f dnA _ -4n,Z,exp(-yJ/(l/y,) 

[dt)~ n m y k z k G(RZ k ) 9 


G(RZ k ) = 



/(I /y k ) = 1 + 0.5./— expOOtl - erf^ k ] 


2 y k 


1 

"4^ 


3 

'8 yt 


where the subscript ^denotes either electrons or ions, 
Tt^is the particle density, Z^= i for the electrons and 1/2 
for ions, R is the magnetic mirror ratio, T^is the charac- 
teristic collision time for electrons or ions, andy^= 
;^e<t>k/7\ where ^is the charge number (= i for singly 
charged ions, = -i for electrons), is the electronic 
charge and 0k is the potential barrier for electrons or 
ions (here and 0e are, taken to be positive quantities). 


The energy loss rates along the field are given by 


d(n k T k ) 

dt 


= z k a/> k 


1 

«My k ) 



A similar estimate of loss rates was made by Sizonenko 
and Stepanov 6 / These equations give the rate at which 
particles or energy are lost by electron or ion diffusion in 
velocity space over their potential barriers <t>e and Oi. In 
view of the expCy^ factor in the particle loss equation, 
the confinement time should increase roughly exponen- 
tially with Computer simulations by McHarg and 
Oakes 68 ^9 show such an exponential increase of parallel 
confinement time with Ok/T^. In the limit of high 
plasma density or low temperature, collisions may be so 
frequent that these ‘Pastukhov’ equations lose validity, 
and the equations of Rognlien and Cutter 7° should be 
used. In MEPC devices the dominant electron losses are 
by diffusion across the magnetic field, but ion confine- 
ment is almost completely electrostatic. 


Parks and Sleeper 1 derived an equation predicting ion 
end loss rates in the transition region between the non- 
adiabatic region (central plasma) and the adiabatic re- 
gion (boundary layer). In the non-adiabatic limit, their 
equation agrees with the result of Tiishmanov7 2 and in 
the adiabatic limit, their equation is similar to the re- 
sults of Chernin and Rosenbluth,73 but with loss rates a 
factor of 3-4 higher, due to effects of magnetic field cur- 
vature. 


2.4 Self-Shielding 


The uncompensated electronic space charge tends to 
shield out the anode potential by an amount A0. As- 
suming azimuthal symmetry, the electrostatic potential 
distribution can be calculated from a two-dimensional 
Poisson equation in the anode gap 


If a_v 

r ydr X dr ) dz 2 


{n-n^e 

£0 


(13) 


Due to the potential hill Oi, only a few ions are passing 
through the anode gaps, so « «l,, as will be discussed 
later. With comb probe measurements in some experi- 
ments it was found that the distribution of electrons in 
the anode regions had bell-shaped distributions across 
the magnetic field. In order to model a variety of condi- 
tions, a Lorentzian-shaped density profile 


(1 + z l !w l ) 


(x 4 ) 


is assumed, where ru a is the electron density at centre of 
the gap and the parameter mo characterizes the half- 
width of the distribution. Assuming this electron den- 
sity distribution in the anode region, with n^ a = 1.5 x 10 ’7 
m"3, a = 2 mm and mo = 0.2 mm, results in the electro- 
static potential distribution of Figure 10. The depth of 
the saddle point relative to the anodes is called A0. If 
the anode length along the magnetic field h » a, then a 
one-dimensional equation 


d<f) (n-n,> , 

— r = (15) 

dz~ £ 0 


can provide an accurate calculation of the potential dis- 
tribution 0(z) across the magnetic field. Assuming 0(-a) = 
0(a) = o, solution of this equation yields 



(16) 


Since high-energy electrons flying through the anode 
gap have trajectories extending one or two Larmor radii 
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on each side of the center, it is assumed that w = 2 p^. 
The resulting values of A0 are shown in Figure 12 as a 
function of n^ a and B, assuming 7T = 20 keV and a = 2 
mm. If ru a = io'9 m 3 and B = 6T, then A 0 = 74 kV. 



Figure 12 — Variation of self-shielding potential sag with magnetic induction, for various 
anode gap peak electron densities n Q , assuming w = 2p Q , a = 2 mm, 7= 20 keV. 


These could be typical parameters for a MEPC reactor. 
A comparable result was found by Ware and Faulkner, 74 
assuming a triangular density distribution. The quantity 
in brackets varies from 0.96 at mula = 0.01 to 0.67 at mala 
= 0.2, so the approximation 




nn a ewa 

2e 0 


(17) 


2.5 Electron Density in 
Anode Gaps 

Five phenomena affect the ratio ions in the gap, 
magnetic reflection, electrostatic acceleration, cold 
trapped electrons and diocotron oscillations. 


2.5.1 Ions in the Gap 

Assuming a quasi-cylindrical plasma with radius r p and 
length L, the rate at which ions are lost is 


- = N2nR28n ia u ia (20) 


where n, is the central plasma density, x, is the ion con- 
finement time, N is the number of ring cusps, R is the 
ring cusp radius, 6 is the effective half-width of cusp 
gaps for ions, n ^ is the density of ions flowing over the 
barrier and is the speed at which ions pass over the 
barrier. Ions just barely leaking over the potential barrier 
will have low velocities m a ~ 3 x io4 m s" 1 , and they will 
be restricted to a narrow corridor of half-width 5 ~ mo. 
At equilibrium the ion particle confinement time will 
equal the electron confinement time, which is roughly 
equal to the diffusion time Hiff across the magnetic field. 
Thus, the ion density in the anode gaps is estimated to 
be 


n ia 

n 


r 2 p L 1 

4RwU ia T diff 


(21) 


gives a slight overestimate of A(f>. Expressing in keV 
and A0 in kV, this becomes 

A</> = 6x10- 15 ^^^(i 8) 

which yields A<jf> = 89 kV (20% high) for the above ex- 
ample case. The maximum electric field occurs at the 
anode wall 



For typical values of MEPC experiments, the electric 
field energy density (i/2)Eoiw* 2 « B 2 /2fi 0 . 


where L, = distance between ring cusps. For the example 
MEPC reactor parameters of table 1, assuming L x = 1.5 
m, it is estimated that n^Jn, ~ 0.02. 

2.5.2 Magnetic Reflection 

Sidorkin and Lavrent’ev75 studied electron reflection by 
multipole cusp gaps, as illustrated in Figure 13. They 
showed that electrons with a a > cw, are reflected by 
the cusp field, where 



x = r/R, r is the initial electron radius, R is the cusp gap 
radius and m is the number of pairs of multipole con- 
ductors = JV/2. The values of cw, are plotted as functions 
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of initial radii in Figure 14, for various values of nu. In 
multipoles with large i the electrons have smaller critical 
angles, so more of them are magnetically reflected. 





Table 1 — MEPC Reactor Parameters 
Example 

The plasma parameters were estimated using approxi- 
mate scalings of attainable density, temperature, con- 
finement time (section 3.7.1) and ^(section 4.1). 

Ring cusp magnetic induction, B a = 6 T 

Applied voltage, (f)A = 400 kV 

Plasma radius. r p = 3 m 

Ring cusp anode radius, R = 4 m 

Ring cusp anode gap half-width, a = 2 mm 


Estimated Parameters 

Equations 

Central plasma density, 
= IO 20 m "3 

(57) and 
Jupiter2M 

Electron and ion temperatures, X^, 
Ti ~ 20 keV 

(48), (49), (58) 

Magnetic field at plasma boundary, 
B b ~~ 1.3 T 

after (37) 

Electron Larmor radius in anode 
gap ,p a = 8 x io _ 5 m 

after (37) 

Electron Larmor radius at the 
plasma boundary , pb = 3.8 x io~+ m 

after (37) 

Self-shielding voltage drop, 
A<|> = 80 kW 

after (47) 

Anode gap peak electron density, 
TL> a = I.3 X IO r 9 m 3 

after (16) 

Energy confinement time, x e j_ = 6 s 

after (40) 

Power gain ratio, Q = 10 

after (60) 


Electrons starting off at rl R = 0.5 near an octupole cusp 
iffL, = 4) have x„ = 0.938 and cw, = 29°. Experimentally 
measured electron beam transmission is consistent with 
such theoretical predictions. Assuming that the central 
plasma electrons have isotropic distributions, the total 
loss fraction through N cusps from a given location is 
approximately 

total loss fraction ~ 0 . 5 ^ [l - cos(a cri( ). ] (24) 

If n, p = density of electrons that are able to penetrate 
into the magnetic cusps, the ratio rupln, can be esti- 
mated from the volume-averaged value of this total loss 
fraction, but this calculation has not yet been per- 
formed. 


Plasma Physics and Controlled Fusion 36 (1994), pp. 1539-1593, doi:io. 1088/0741-3335/36/10/001 


H 



2.5.3 Electrostatic Acceleration 2.5.5 Diocotron oscillations 


Electrons that succeed in entering the cusp gaps are ac- 
celerated by the positive applied voltage there. Moir et, 
ah 6 showed that the 



electron density in the gaps is lower than the magneti- 
cally penetrating density ru p by a factor 


— = exp {y)erfc(yjy) = 

h p 



3 15 

jny 




(25) 


where y = 0i/T t ,, and erfc is the complementary error 
function.77 For values of = 3 to 6, this ratio is n^Jn^p 
= 0.29 to 0.22. 


2.5.4 Cold Trapped Electrons 

Magnetic reflection and electrostatic acceleration tend 
to reduce but electrons produced by ionization of 
neutral gas tend to become electrostatically trapped in 
the anode region (orbit 4, Figure 11), increasing If 
classical cross-field diffusion were the only loss process, 
then the accumulation of trapped electrons could in- 
crease n^a and A </> to unacceptably large values, but the 
diocotron instability may help remove cold trapped 
electrons without seriously impairing hot electron con- 
finement. 


The electric field due to electron space charge is zero at 
the middle of a ring cusp gap, increasing on each side. 
The resultant ExB drift velocity varies spatially, and the 
shear of this velocity gives rise to the diocotron (slip- 
stream) instability" 8 ? 1 ; Computer simulations by Fevy 
and Hockney 80 showed good agreement with theoretical 
instability growth rates. Pankrat’ev etal Sl found that the 
long-wavelength instability (X> io'w, where mo is the elec- 
tron stream half-width) occurs when the electron drift 
velocity is near the wave phase velocity, and it is sup- 
pressed by a nearby conducting wall. However, Gordi- 
enko etal 81 found that the presence of two groups of 
electrons (trapped and passing) in the anode gaps broad- 
ens the unstable range of wavelengths and makes dioco- 
tron oscillations unstable whenever there is a vacuum 
gap between the electron stream and the anode wall. 

The short-wavelength instability (k < 0.2 a) occurs near 
the electron cyclotron frequency, 



and has a growth rate where q = c w pi J-l "WceX- These dio- 
cotron oscillations play a beneficial role by removing 
cold trapped electrons from the anode region, where 
they would broaden the width mo of the electron density 
distribution and increase A </>. The instability growth rate 
is small if q < 0.2, which is equivalent to the condition 
that 


n a <(2jtl0 18 m _3 )B a 2 ) 

where B a is the magnetic induction in the gap (T). Thus, 
as TL, a decreases well below this limiting value, the short- 
wavelength diocotron instability turns off! At B a = 5 T, 
this limit is < 5 x xo : 9 m'3, which is less restrictive 
than the value of il* < 1.4 x ro'9 m'3 required to keep A </> 
< 100 kV (Figure 12). 


Yushmanov 8 3 studied the build-up and removal of 
trapped electrons in the anode regions, and determined 
the required anode dimensions. The width h of the an- 
ode (Figure 11) should be in the range 


5asisi - 56 ' , (®)(S 8 


where a - 4.7, e = trapped electron pump-out rate/ 
diffusion rate = 1/3, N is the number of cusp gaps and 
rujrua is the trapped electron density/anode electron 
density = 1/3. Taking N = 6 in the example case of table 1, 
it is found that 
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h m ax = 4 cm. Low-frequency oscillations can selectively 
pump out the trapped electrons without seriously in- 
creasing the loss rates of the free plasma electrons. 

If the electrons in the central plasma have a quasi- 
Maxwellian distribution, the distribution will be up- 
shifted in energy by an amount e0; as the electrons enter 
the anode regions, as illustrated in Figure 15. Electrons 
produced by ionization in the plasma edge region will 
have a spread of energies. The diffusion rate is high for 
low-energy electrons, due to interaction with diocotron 
oscillations, and it may be nearly classical (D 0= WeJ l/z ) 
for high-energy electrons. The resulting electron energy 
distribution in the anode regions is expected to resem- 
ble the shifted Maxwellian distribution of Figure 15 
(middle). 



distribution in anode region, assuming 4> / = 2T g (middle); relative diffusion rate in anode 
region, with effect of diocotcon oscillations on cold electrons (bottom). 

Some magnetically trapped electrons have low enough 
energies that they can immediately interact with dioco- 
tron oscillations; but most of them, having been acceler- 
ated by the anodes, must be decelerated by collisions 
before the diocotron oscillations interact strongly with 
them. For one example case s 4 the trapped electron frac- 
tion was estimated to be rujn^a ~ 0.24. Considering ions, 
magnetic reflection, acceleration, cold trapped electrons 


and diocotron oscillations, there will probably be about 
an order of magnitude density ratio between the anode 
region and the central plasma. 

^~ 0.1 

n (29) 

The ratio would probably be lowest for cases with a 
large volume of field-free plasma, narrow anode gaps and 
low neutral gas pressure; and it could be near unity for 
cases with spindle cusp magnetic field, wide anode gaps 
and high neutral gas pressure. Values inferred from data 
in a variety of experimental conditions range from 0.01 
to 1. 

2.6 Plasma Potential 

At Te, > x keV and ru < xo 20 mf recombination is negligi- 
ble. Ions are confined electrostatically with energy loss 
by charge exchange and by diffusion in velocity space 
over the potential barrier. The ion and electron particle 
conservation equations may be written: 

dn I _ ^ ( dn t \ 

dt V dt J< 30 ) 

dn _ S + I (dn\ n 
dt i + eV { dt J,| t x(3i) 

where S, is the ionization source, I is the electron cur- 
rent from cathodes into the plasma and Xj_ is the elec- 
tron loss time by diffusion across the magnetic field. At 
equilibrium the plasma potential will adjust itself so that 
these electron and ion loss rates are equal 8 ? 



Since 

<l>e =<t>A~ 4 >i ~ A( t> ( 33 ) 

these two equations can be solved simultaneously 86 for 
the two unknowns and (j>i using equation (9) for the 
parallel loss rates. The plasma potential distribution has 
been computed self-consistently in some two- 
dimensional models. 

2.7 Two-Dimensional Models 

Brunei etafii solved a two-dimensional (r, ^ Poisson 
equation, assuming nearly Maxwellian distribution func- 
tions dependent on the constants of motion. They cal- 
culated the distribution of potential and particle densi- 
ties in an electrostatically plugged spindle cusp corre- 
sponding to the KEMP device. They found that the 
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electron and ion densities were both maximum on-axis. 
The electron density had a secondary peak in the ring 
cusp anode, and the ion density was very low in the an- 
ode regions. Matte and Lafrance 88 developed an im- 
proved numerical method for solution of the Poisson 
equation in MEPC devices. 

Shoucri etafi 9 simulated an MEPC plasma with about 
5000 electrons and ions (M/zzl = 16, where M is the ion 
mass and m, the electron mass) in a 128 x 50 cylindrical 
grid, using a capacitive matrix Poisson solver, with pa- 
rameters similar to those of the KEMP II experiment. 
The simulation showed the maintenance of a two- 
dimensional electrostatic potential well, with saddle 
points in the anode gaps, as expected. Electrons were 
very well confined, but ion confinement was limited by 
plasma self-shielding in the point cusp gaps, especially 
during electron injection. The confinement time ob- 
served in the simulation was scaled up by the square 
root of the experiment/simulation density ratio to yield 
a predicted confinement time of 50-75 ps for the KEMP 
II experiment, which was consistent with experimental 
measurements. 


2.8 Energy Confinement Time 

The non-radiative cross-field electron energy loss time 
may be written 

_L-_L 1 


where w is the time for heat loss by conduction -1.51 
t diff, Tdiff is the characteristic time for diffusion across the 
magnetic field and t trap is the time for free electrons in 
the central plasma to become trapped by the boundary 
magnetic field. 


One estimate of the trapping time is? 0 


2 mvB t 
V rce 3 In A 


( 35 ) 


where nj is the electron velocity and In A is the Coulomb 
logarithm. For reactor conditions T trap « ~X diff, and 


T e± -0.6T diy 


(36) 


2 T ei aV 

S0C PbPa (J7) 


where T„ is the electron-ion momentum-transfer colli- 
sion time, V is the plasma volume, S is the plasma sur- 
face area, pb = (2mkTe) ll2 /eBt is the electron Larmor ra- 
dius at plasma boundary with B/, = {2p 0 nk^X^ + T,) l/2 , p a = 
(zmkTeJ I,2 leB a is the electron Larmor radius in the anode 
gap and a = a parameter varying from 1.33 at 
T , 7 Te, = o to 4.6 at TilTe. = x. A similar estimate was de- 
rived by Dolan etal. 9 2 Assuming In A = 16 in a hydro- 
genic plasma, the electron-ion momentum-transfer col- 
lision time 44.93 may be reduced to the form: 

_ 9 . 4 x 10 14 r g f 


where T^is in keV For spindle cusps, V/SR ~ 0.01-0.03, 
depending upon dimensions, plasma heating and con- 
finement effectiveness. (Low-pressure plasma have 
smaller volumes.) This low ratio plus diffusion enhance- 
ment by electron collisions with neutral gas give spindle 
cusps poor confinement times. 


Assuming classical diffusion in a spindle cusp with radius 
R and length L, Lavrent’ev^ derived analytical expres- 
sions for the electron density profile, the plasma volume, 
the average electron density and the electron confine- 
ment time. The electron energy loss rate is: 


P e = [ 0 ,. + A(f) + 2 T e (cosh (p) - 1 )] 


nR 4 D 0 n 0 p 

r a 2 Lsinh(p) 


(39) 


where D 0 is the collisional (classical) diffusion coeffi- 
cient, is the central plasma density, p = (x±/xt||) i/2 . and 
r a is the point cusp anode radius. This power loss is 
minimized when p - 5. Lavrent’ev found that the edge 
electric field can cause significant increases of electron 
transport rates in. small devices with substantial neutral 
gas pressures. For the parameters of the JupitenM de- 
vice the electric field and electron mobility should result 
in about a 45% increase in the electron cross-field trans- 
port above the classical diffusion rate. In the case of a 
reactor with ‘burn-out’ of the neutral gas, however, the 
electron mobility would be small, and the electric field 
effect would be diminished. 


Nearly classical transport has been observed in some 
MEPE experiments, while rapid anomalous transport 
has been seen in others. 

2.8.1 Classical Transport 

Pastukhov? 1 estimated the classical diffusion time to be: 


Lavrent’ev^ showed that angular momentum from ExB 
rotation produces an effective potential well ffy that 
greatly inhibits axial ion losses, even if the point cusp 
electrostatic barrier = o. A similar effect acts on elec- 
trons, but is rn^lM times weaker. Near the ring cusp the 
electric field effects dominate the ion motion. Dol- 
gopolov etal s’ 6 computed ion trajectories in an electro- 
statically plugged spindle cusp, and found that the mag- 
netic field gradient force and the centrifugal force (re~ 


Plasma Physics and Controlled Fusion 36 (1994), pp. 1539*1593, doi:io. 1088/0741-3335/36/10/001 


17 


lated to the azimuthal ExB drift) tend to keep ions 
from escaping through the point cusps. This inhibition 
of point cusp ion losses makes plasma confinement ef- 
fective in a linear set of ring cusps. 

For a linear set of ring cusps, the volume-to-surface area 
ratio VISR = o. 1-0.3 * s 111 order of magnitude better than 
or spindle cusps. Assuming VISR = 0.02 and 7V, = T„ 
equation (37) becomes 


(40) 


For the example MEPC reactor parameters case (table 1) 
it is found that x„ = 0.84 ms, the diffusion time Xjijf = 20 s 
and x el = 12 s. Since experiments have attained about 
half the classical value, a value of Xel = 6 s is fisted in ta- 
ble x. This near classical transport would make an attrac- 
tive reactor, but anomalous transport could prevent a 
reactor from being successful. 


2.8.1 Anomalous Transport 

Pastukhov and fl’gisonis97„98 studied anomalous electron 
transport in the boundary layer of an MEPC plasma. 
When Ti « 7V. , the ion acoustic instability is dominant, 
giving rise to rapid electron diffusion through the 
plasma edge layer. The long-wavelength ion-acoustic 
instability can occur at k - 1 and w - (tw C£ ,co„) 0 -5 when 
the electron density-gradient drift velocity ■x'vn > c s . The 
corresponding maximum diffusion coefficient D max - cp^, 
which has Bohm-like scaling. The resulting anomalous 
transport tends to broaden the electron density profile 
to a condition where 'X'vn - c s . The resulting theoretical 
electron density profile is given by: 

n(r ) = n 0 exp (~ec s A / cT e ) ( } 

where ru 0 is the central plasma density and A is the 
magnitude of the magnetic vector potential. The waves 
are partially absorbed by the ions, and may heat them to 
the point where the instability turns off. When T ,■ - 7V,, 
the ion acoustic instability becomes a Buneman instabil- 
ity, and the instability growth rate decreases by a factor 
of about (pz,/M) l/2 . Once electrons diffuse outside the 
bounding flux surface (the surface which just touches 
the anodes), they are magnetically mirrorconfined (orbit 
3 of Figure 11) until they scatter into the loss cone and 
travel along the magnetic field to the anodes. For a fin- 
ear quadrupole with alp^< 2 (MlmJ) lh , the anomalous 
electron confinement time is 



where X 0 is the electron collision time. This case is simi- 
lar to that of the ATOLL experiment, and the resulting 
confinement time is much shorter than with classical 
diffusion. Anomalous transport would spoil the energy 
balance in a reactor, resulting in a low fusion power den- 
sity 


2.9 Energy Balance 

As neutral atoms impinge on the plasma some are ion- 
ized in the boundary layer. The resultant ions are accel- 
erated as they fall down the potential hill into the 
plasma region. However, Pastukhov9 I found that this ion 
heating process is inefficient, and that ion heating by 
electrons is very slow. Ions are also heated by plasma 
waves. Ion energy is lost by charge exchange and by con- 
vection, as heated ions escape over the potential hill 
Since ions with energies Wi > eO, are lost out the cusps, 
their distribution is truncated at higher energies, and 
the resulting fusion reaction rate is reduced.99 

Electron energy may be supplied by electron beam injec- 
tion, neutral beam injection, wave heating and heating 
in the edge electric field. Electron energy is lost by elec- 
tron flow over the potential barrier to the cathodes, by 
transfer to waves, by thermal conduction and convection 
across the magnetic field and by inelastic collisions with 
ions, atoms and molecules. 


Approximate energy balance equations for the central 
plasma electrons and ions may be written 


L5 ^ = ^ + rdWl 

dt eV E ' X ' a ae L dt 

-p 


-1.5 — --P.-l.5n— 


(43) 


(r e -7;) 


+ f l S m W a + F 


+0.25 nf (< 3 v)(f ai W a - 1 .8e</> ( ) 

v dt A\ (44) 


where Pe is that portion of electron heating by the edge 
electric field which affects the central plasma, x eq is the 
electron-ion equipartition time, P n d is the radiative 
power loss, Wi, m is the average electron energy loss per 
incident neutral atom due to ionization, P ^ is the elec- 
tron power loss to plasma waves, f is the fraction of in- 
cident neutrals which are ionized (instead of causing 
charge exchange), f a ^ and fai are the fractions of fusion 
product alpha energy W a which are transferred to elec- 
trons and ions, S m , is the volume-averaged neutral input 
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rate, W a is the average energy acquired along the slope 
of the potential well by ions, <O r v> is the DT fusion re- 
action rate parameter and Pi is the ion heating by plasma 
waves (such as ion acoustic). It is assumed that the aver- 
age energy of each fuel ion undergoing a fusion reaction 
is about o.yeJbi. Since the alpha particles are poorly 
confined, fa^« i and f a i « i. Additional terms would be 
needed for neutral beam injection, microwave heating, 
etc. These equations are nonlinear, because several pa- 
rameters depend on the density and temperatures. Simi- 
lar equations have been solved numerically for several 
cases. 8 t Typical conditions satisfying the plasma poten- 
tial and energy balance equations were found to be 


t 

~ 0.5 


<t>A 


(45) 

01 

- 0.3 


0a 


(46) 

A</> 

- 0.2 


0A 


(47) 

T^0.05(j) A 

(48) 

T e ~ C 

).05(j) A 

(49) 


In order to confine plasmas with fusion reactor parame- 
ters, applied voltages 0a 2 300 kV are desirable. 100 


core is essentially spatially isotropic, but the electrons 
do not have a Maxwellian velocity distribution. The ion 
flow is expected to remain spatially anisotropic in the 
mantic region and to retain its non-equilibrium (nearly 
mono-energetic) nature over typical ion lifetimes. 


Startup 


Startup Beam Injector 
Through Point Cusps 





Dense Plasma 
Core with 
Radius r 


High-Beta Plasma Operation 


Karpukhin etal 101 studied particle and energy balance in 
an electromagnetic trap with r f - 0.3 m, B a = 1-3 T, 0a = 
25-50 kV, V - 0.5 m 3 , a = 1.5 mm and an electron injection 
power of 500 kWi They predicted values of ru - 3 x xo'9 
m _ 3, te - 0.1-0.2 s and T; - 0.5-1 keV. For an octupole elec- 
tromagnetic trap they found the optimum electrode gap 
spacing a = 0.5-1 mm, which maximizes the injection 
efficiency, plasma density and confinement time. 102 En- 
ergy balance calculations will be discussed further in 
section 4.1. In addition to spindle cusps, multipoles and 
ring cusps, two additional MEPC configurations have 
been studied recently: the Polywell concept and the 
modified Penning trap. 

2.10 The Polywell Concept 

Bussard to 3 proposed electrostatic confinement inside a 
polyhedral array of magnetic point cusps (Figure 16), 
called a Polywell, with electron injection at 10-100 kV 
along point cusp axes to provide a quasi-spherical nega- 
tive potential well capable of trapping and confining 
ions. The ions are to be injected at very low energies (< 
100 eV) from the system boundary (r = R ), along cusp 
axes or the axes of the polyhedron vertices, The electron 
distribution in the ‘mantle’ region (r > r c ) outside the 


Figure 1 6 - Electron injection through multiple point cusps in a Polywell device. The full 
curves represent magnetic field lines, which are pushed outward by electron diamagne- 
tism during high-beta operation. 109 

The degree of convergence of ions flowing radially to 
the central core is determined by angular momentum 
conservation to be 



where Wl is the ion energy perpendicular to the radial 
direction and 0 i is the accelerating potential . I0 4 The 
Polywell concept is based on high-beta electron con- 
finement by the magnetic field. Watrous etal 10 ^ showed 
that the required diamagnetic currents could be 
achieved with reasonable electron injection conditions, 
yielding a relatively field-free internal region that en- 
hances electron confinement. Krall etal 106 also studied 
the density and pressure anisotropy profiles resulting 
from various ion sources, and assessed processes which 
may degrade ion focusing, such as angular deflection and 
energy upscattering. The concluded that isotropy of the 
plasma core inhibits the Weibel instability 
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Solving the Poisson equation for spherically convergent 
charged particle flows, King and Bussard I0 7 found that, 
with a high degree of electron recirculation within the 
system (with long electron lifetimes), central plasma 
densities of about ro 2 ' m"3 could be attained at high elec- 
tron gun currents. Bussard and King 108 distinguished 
two bounding operating regimes: a low-electron-beta 
mirror regime, where magnetic trapping and reflection 
are dominant, and a ‘whiffle-ball (WB) regime’, where a 
high-electron-beta plasma pushes the magnetic field 
outwards, so that the plasma is situated inside a quasi- 
spherical magnetic surface with small holes correspond- 
ing to the point cusps. In the WB regime they find that 
the electron recirculation factor Gj 0 (representative of 
the central plasma electron lifetime) increases almost 
linearly with /L c r c 2 , where ru c is the central electron and 
ion density within the core radius r c . They advocated 
plasma start-up at low magnetic field, in order to avoid 
the mirror regime, which requires higher electron beam 
currents. 

Three issues arise with the Polywell concept: 

(1) It may be difficult to achieve a high degree of 
ion focusing (small W±), due to ion reflection by 
the scalloped electrostatic potential well 
boundary I0 9 and ion deflection by wave-particle 
interactions. Larger W± would result in a larger 
core radius r c and reduced central plasma den- 
sity 7L C . Further research is needed to clarify the 
attainable degree of focusing. Even if great fo- 
cusing were not achieved, the Polywell might 
still function as a fusion reactor by using a 
larger radius R. 

(2) Initially mono-energetic electron or ion distri- 
butions may be driven toward thermodynamic 
equilibrium (approaching Maxwellian distribu- 
tions) by enhanced collisional mechanisms, such 
as the beam-plasma and two-stream 
instabilities. 00 It would be desirable to prevent 
the accumulation of cold electrons, in order to 
avoid the beam-plasma instability. 

(3) The polyhedral set of point cusps has magnetic 
field nulls on its surface, which might not con- 
fine high-pressure plasma well. However, no 
significant particle losses were associated with 
the null points in the early operation of the 
spherically convergent ion focus (SCIF) 
experiment. 111 

2.11 Modified Penning Trap 

Barnes and Turner 112 found that a time-varying external 
electrostatic quadrupole field in a uniform static mag- 
netic field could induce large-amplitude, nonlinear oscil- 
lations of a non-neutral plasma with very high plasma 


densities at the centre in one-, two- and three- 
dimensional cases. These oscillations are predicted to be 
stable at large amplitude. 

By definition, the Brillouin ratio is the ratio of relativis- 
tic electron energy density to magnetic field energy den- 
sity: 

„ . . 2 u n nmc 2 n 

Brillouin ratio = 3 = — 

B- n, (;l) 

where tub is the Brillouin density. In an ordinary equilib- 
rium case £ 1, the ‘Brillouin limit’. Turner and Bar- 
nes °3 found that a class of flows with deviatoric strain 
can support charged particle densities with Brillouin 
ratios locally exceeding unity This configuration does 
not require large-amplitude plasma oscillations. They 
showed that, in a toroidal geometry with a purely poloi- 
dal magnetic field, local values of the Brillouin ratio n^l 
tub>i could be sustained. The poloidal magnetic field 
could be produced by a toroidal quadrupolar current 
configuration outside the plasma, and a toroidal octupo- 
lar electric field applied outside the plasma could be 
used to sustain the required plasma configuration. With 
ai/r term in the toroidal flow velocity, the Brillouin ra- 
tio at the x point (the center of the plasma) becomes 
infinite. In contrast to the usual Penning-trap-confined 
non-neutral plasmas, this configuration is not in ther- 
modynamic equilibrium. Large values of the Brillouin 
ratio are feasible in a cold plasma when the electron ve- 
locity shear length L, ~ •v/wp^. 

Barnes etal IJ 4 studied non-neutral plasma confinement 
in a uniform magnetic field with an applied electrostatic 
quadrupole field (a Penning trap). They found a class of 
confinement that does not require a large velocity shear. 
The plasma reactivity is proportional to <n 2 >, where < > 
denotes a space and time average. Since the Brillouin 
limit applies only in the volume-averaged sense, the re- 
activity can be enhanced by inducing large spatial varia- 
tions of ru over the plasma volume (such as ultrahigh hl, 
at the centre). Low-energy, low-Pe (canonical angular 
momentum) charged particles (such as electrons or deu- 
terons) could be injected into this spherical Penning 
trap. The effective parabolic spherical potential well of 
the applied electromagnetic field would accelerate the 
particles radially inwards to energies -100 keV. Conver- 
gence near the centre (r = o) produces a dense, inertially 
confined, non-neutral plasma core. Strong focusing oc- 
curs only for a single charge-to-mass ratio qlnu, so a pure 
deuterium plasma might be used. Miyamoto el al °5 
showed that, with a judicious choice of a rational integer 
determining the electric field shape, the DD fusion re- 
action products 3 H and 3 He could also be focused to 
pass through the central region, enhancing the fusion 
power density 


Plasma Physics and Controlled Fusion 36 (1994), pp. 1539-1593, doi:io. 1088/0741-3335/36/10/001 


20 


Another possibility is an initially pure electron plasma, 
with deuterium-tritium ions trapped in the virtual cath- 
ode. In the electron-plasma version low energy electrons 
would be introduced through one of the negative end- 
caps of the electrostatic quadrupole, as illustrated in 
Figure 17. In a reference frame rotating at angular fre- 
quency Q = -eBlzm,, the Poisson equation becomes 1 ^ 

V 2 0 = f-L B -n) 

V e J (52) 

The applied electric field is produced by electrodes with 
hyperbolic surfaces tangent to the spherical boundary of 
the potential well. 



Particles with energies <zV 0 l^ will be confined, where 
the applied voltage V 0 = ruRea 2 ! \z a . 

The cold injected electron beam will be reflected near 
the centre of the spherical well (r = o). If Pe = o, then the 
characteristic radius r c of the dense central core could be 
comparable to a local Debye length. With finite perpen- 
dicular injection temperature T_ the angular momentum 
conservation equation yields 



where W is the electron injection energy and r 0 ~ a !6 is 
the radius of a virtual anode potential peak in the rotat- 
ing reference frame. The electron beam velocity nj varies 
slowly, except near the inner and outer turning points. 
From the steady-state continuity equation with Pe = o, 


the density varies radially as nJr) = nXd)a 2 lr 2 , where a is 
the radius of spherical bounding surface, and the central 
core density ru c = ruB&/r c 2 . Since nJ 2 decreases rapidly for 
r > r c , the average of nJ 2 is estimated to be 


c (54) 


This proportionality to V„lair c shows that the reactivity 
can be enhanced by high voltages and small radii. 

The scaling was confirmed by numerically computing 
self-consistent Vlasov equilibria for an example case 
with a = 3 mm, V 0 = 120 kV, electron beam radial injec- 
tion at W = 50 eV, I = 0.125 A, Ti_ = io _ 4 eV and T r = 5 eV. 
The numerical result showed that a density = 4 x 10 2 J 
m’3 could be confined within r c ~ 0.4 um.°4 For compari- 
son, the above approximate equations yield r c = 0.71 um, 
rLjB ~ 3 x 10 18 m _ 3, n^c ~ n^na 2 lr c 2 = 5.3 x io 2 ? mf3 and <tl, 2 > = 

3.7 x io4 m 3. Thus, for this example case the approxi- 
mate scalings yield a core radius and density not far from 
the numerically computed values. 

A second example case considered the effects of ions 
that partially neutralize the electron space charge. As- 
suming tzw/tl = 0.5 for r < 1 um, hw/hl, = (1-10 6 ) for 1 um < 
r < r 0 , ruiln, = o for r>r 0 and I = 0.425 A, similar poten- 
tial and density profiles were obtained, with r c = 0.3 mm 
and = 4 x io 2 5 m"3. The neutral background penetrates 
most of the plasma, but within a small ‘burn-out’ radius 
high ion temperatures may be sustained. The central ion 
confinement time is about 0.25 ms, and a neutral gas 
pressure about 1-10 mPa would be appropriate. For this 
example case the predicted central ion temperature was 

7.8 keV, and the DT fusion power would be about 50 
mW, and the energy gain Q = io"3. Calculations indicate 
that this configuration is stable against the two-stream 
instability, that Coulomb collisions do not significantly 
alter the electron beam velocity distribution, and that 
small angular momenta are attainable. :i 4 

Although such a Penning trap would not produce a large 
total power output, it could produce useful fusion reac- 
tion rates in a laboratory scale device. Higher-order mul- 
tipole traps could generate higher power outputs. 11 ^ 

Tiouririne etal 21 1 studied the use of higherorder multi- 
pole traps for non-neutral plasma confinement. The to- 
tal electrostatic potential 0 is the sum of the externally 
applied potential 0 a pp and the potential 0 S generated by 
the plasma space charge. The plasma density, computed 
in terms of 0 using the Poisson equation, was assumed 
to be zero outside a closed boundary, and the plasma 
source potential 0 s was computed assuming that it van- 
ishes at large distances. Then the required value of ap- 
plied potential 0 a pp = 0 - 0 s. The magnetic field contours, 
plasma density contours and externally applied potential 
contours computed for the magnetic quadrupole Q = 2) 
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and magnetic hexapole (J = 3) cases in cylindrical geome- 
try are shown in Figure 18. 


1 = 2 (magnetic quadrupole) 
(a) isodensity 





l = 2 (magnetic quadrupole) 
(b) magnetic flux 



(c) externally applied 
electrostatic potential 



Figure 18c — Contours of externally applied electrostatic potential for a magnetic quadru- 


1 = 3 (magnetic hexapole) 
(d) isodensity 
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I = 3 (magnetic hexapole) 
(e) magnetic flux 





(f) externally applied 
electrostatic potential 



Figure 18f — Contours of externally applied electrostatic potential for a magnetic hexa- 
pole. 117 

The global Brillouin ratio was found to be ^ = zlKJ + x), 
which approaches z for high-order multipoles. 

Turner etal 118 calculated the frequency of the non- 
neutral plasma anharmonic oscillations and found that 
(olojc varies from z~ lh when 5 = o to (zlj) 1 ' 1 when 8 = 1, 


where 8 = (p, -pi)l(p 2 + pi), and pi, pi are the turning 
points of the oscillations of the plasma boundary 

Although the studies of modified Penning traps are just 
beginning, they have shown several interesting features: 

(i) In a non-neutral plasma, densities exceeding the 
Brillouin limit can be confined by several 
means: in an oscillating plasma; in a flow with 
deviatoric strain (related to velocity shear) such 
as a torodial magnetic quadrupole with applied 
electric field; in a homogeneous magnetic field 
with an applied electrostatic quadrupole and 
low-energy particle injection; in a higherorder 
magnetic multipole field with an applied elec- 
trostatic multipole. 

(ii) A large degree of focusing is possible, with re- 
sultant very high densities inside the core 
radius tv. This focusing can result in high values 
of reactivity, which is proportional to <nJ>. 

(iii) Higherorder magnetic multipoles can have 
global Brillouin ratios approaching 2. 

While these results are not immediately conducive to an 
attractive fusion power reactor, further developments 
along these fines could result in newer configurations 
with more attractive properties. The evolution towards 
higher-order magnetic multipoles is similar to that of 
the electromagnetic trap experiments, which started 
with spindle cusps, then went to the toroidal quadrupole 
and higherorder multipoles, such as a linear set of ring 
cusps. 
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3. MEPC 
Experiments 

3.1 Spindle Cusps and Mirrors 

Early experimental studies of electrostatic plugging were 
conducted by Lavrent’ev*’* 6 and by Hilton etal. uc > Due to 
financial limitations, most early experiments were 
pulsed spindle cusps, although spindle cusps are poor 
configurations for electrostatic plugging. Electrons were 
injected from electron guns in the point cusps, with 
plasma production by ionization of neutral gas. Pulsed 
magnetic fields have several disadvantages for MEPC 
experiments: 

(i) Eddy currents cause magnetic field errors; 

(ii) During magnetic field decay the plasma drift is 
outwards; 

(iii) Plasma density build-up and heating may take 
longer than the magnetic field lasts. 120 

Experiments with spindle cusp electromagnetic traps 
like S-i, S-3M, S-4M, and JupitenA showed that: 

(i) The electron confinement time increased by 
three orders of magnitude when the plugging 
voltage was applied. 

(ii) The electron confinement time was inversely 
proportional to the neutral gas pressure. 121 

(iii) The rate of electron loss across the magnetic 
field was similar to that expected from classi- 
cal diffusion due to collisions with ions and 
neutral atoms. 

(iv) The plasma density increased linearly with 
applied voltage, up to a stability limit. 

(v) The plasma density increased proportional to 

(vi) The mean ion energy scaled linearly with ap- 
plied voltage . I2 4 

(vii) The electron distribution function 81 shown in 
Figure 19 is comparable to the theoretical 
curve of Figure 15. 



W, keV 


Figure 19 — Measured electron energy distribution in the anode region immediately after 
cut-off of electron beam injection, and 1 .5 ms later.81 

(viii) A deep, negative potential well could be sus- 
tained for many ms, provided that the neutral 
gas pressure was low. 

(ix) Electron orbits and loss,, angles agree with 
theoretically predicted values, provided that 
the starting point is within an adiabatic 
region. 12 * 

(x) Diocotron oscillations occurred, as predicted 
theoretically 

Experimental observations of diocotron oscillations in 
the ring cusp indicated the following: 

(1) The frequency of the oscillations was propor- 
tional to applied voltage and inversely propor- 
tional to magnetic field strength, as expected. 

(2) The growth rate of the instability was propor- 
tional to neutral gas pressure: and thus to the 
density of cold electrons produced by ioniza- 
tion. 

(3) By segmenting the ring cusp cathode, the rela- 
tive ion currents as a function of azimuthal an- 
gle were measured. These measurements indi- 
cated that the long-wavelength instability had a 
mode number »L = ior two wave periods 
around the circumference of the ring. Higher 
mode numbers were apparently more stable. 126 

(4) The anode region electron density remained 
below the short-wavelength diocotron limit of 
equation (27). 

Sugisaki l2 " showed that the application of electrostatic 
fields to a magnetically confined plasma prolonged the 
plasma decay time. Strijland 128 studied particle orbits in 
an electrostatically plugged spindle cusp and the varia- 
tion of confinement time with diaphragm diameter. 
With electrostatic plugging of a magnetic mirror, 
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Nishida etal 126 observed a reduction of ion losses to 15% 
of the unplugged value, and ion heating by E x B rota- 
tion. 

Vdovin and co-workers '3° studied electron loss processes 
from the S-4M electromagnetic trap. They derived an 
equation for equilibrium plasma density that varied with 
electron beam current, electron energy and neutral gas 
pressure in accord with experimental data. Plasma elec- 
tron diffusion was dominated by scattering off gas mole- 
cules and molecular ions. They suggested that anoma- 
lous electron transport in the ring cusp region may be 
caused by a long-wavelength two-stream instability. The 
injected electrons had an anomalously fast loss rate in 
the S-4M electromagnetic trap, and neutral gas effects 
kept the plasma cold (a few eV temperature). 

According to measurements by Komarov et ah' with a 
two-wavelength microwave interferometer in the 
Jupiteri A spindle cusp device, the central density was 
an order of magnitude higher than the average density. 
Zalesskii et al 1 ? 2 found that the the plasma accumulation 
rate with microwave heating was much faster than with 
electron beam injection, and about 50% of the micro- 
wave power was coupled to the plasma. 

Stepanenko and Komarov 1 ?? measured the Jupiteri A 
electron and ion energy spectra during microwave heat- 
ing with 10% power absorption, and determined the 
electron energy confinement time to be about 0.9 ms. 
The dominant ion energy loss mechanism was charge 
exchange. 

Electron and ion energy spectra were measured with 
gridded electrostatic energy analyzers at the point cusps. 
The electron energy distribution is shown in Figure 19. 
The ion energy spectra were measured with a drift mass 
spectrometer combined with an electrostatic ion-energy 
analyzer. The distribution appeared to have two Max- 
wellian components. The higher temperature compo- 
nent was believed to be caused by the initial formation 
of the potential well. The average energy of the ions 
increased linearly with applied voltage, <W,> - 0.18 <£a.? 6 

Komarov and Stepanenko 1 ? 2 ? measured the spatial varia- 
tion of the total number of electrons AC leaving the ring 
cusp gap, and found a narrow distribution with width 6 - 
p^, as shown in Figure 20. At an electron energy of 0.55 
keV, pe_ = 0.1 mm. The distribution was narrower at high 
plasma densities, indicating that electron losses along 
the magnetic field became important relative to cross- 
field diffusion losses. The model profile with r w =p<^ is 
roughly consistent with the experimental data, but -w = 2 
pt_ is used here in estimating future device parameters to 
avoid underestimating A 0 . Komarov and co-workers‘35 
also measured the distributions of electrons and ions in 
a point cusp gap, finding a narrow electron stream with 
radius <1 mm. 



Figure 20 — Experimentally measured electron density distribution (full curve) in the ring 
cusp of the Jupiter 1 A device at cD^ = 6 kV, B = 0.82 T, and model profiles nln a = 1/[1 + (z / 
w)2] with w = P c (broken curve) and w = 2p c (dotted curve). 

In the KEMP-II spindle cusp experiment Stansfield and 
colleagues 1 ? 6 measured bl, - 4 x io'7 nr?, and electron en- 
ergies - 1 keV. Ionization of gas prolonged the plasma 
density persistence after electron beam cut-off. 

3.2 Jupiter-1 M 

The JupitenM spindle cusp experiment at Kharkov 1 ?? 
had pulsed magnetic point cusp fields < 13 T (25 ms half- 
period), ring cusp fields < 4 T, electron injection voltages 
-4 kV, plugging voltages 0a- 6 kV, electron beam injec- 
tion current <2 mA, neutral gas pressure 0.1-1 mPa and 
plasma volume 30 cm?. It attained plasma densities up to 
10 18 m'?, with temperatures -100 eV and confinement 
times - 1 ms. The short duration of the pulsed magnetic 
fields did not allow enough time for build-up of higher 
plasma densities. The plasma losses from the point 
cusps were less than from the ring cusp. 1 ? 8 

Komarov and colleagues 1 ? 6 used gridded analyzers to 
measure the energy distributions of electrons and ions 
emerging from the cusp gaps after removal of the plug- 
ging voltage. The mean ion energy observed under those 
conditions was about one-fifth of the plugging voltage. 
For a few microseconds after electrostatic plugging was 
removed, the remaining ions were accelerated to higher 
energies by the plasma potential. 12 ? 0 The measured po- 
tential barrier in the ring cusp for ions 0j - 3 T, exceeded 
the calculated barrier - 2 T;. The barrier height in the 
point cusp was practically zero, but most of the ion 
losses occurred through the ring cusp. Point cusp ion 
losses were probably retarded by the centrifugal and 
magnetic field gradient forces. 12 ? 1 
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Azovskii and colleagues 1 * 2 found that the inner bound- 
ing surface of the plasma in the JupitenM experiment 
has a radius r„ = (2 Rp 0 ) 112 in the point cusp (where p 0 is 
the electron Larmor radius in the point cusp), as pre- 
dicted from flux mapping to the ring cusp. The density 
of plasma flowing out of the point cusp was two orders 
of magnitude lower than the central plasma density, and 
A0~ 0.5 kV was inferred from the shape of the ion en- 
ergy distribution emerging from the point cusp .'43 

Azovskii and co-workers 1 ** measured the plasma density 
in Jupiter iM simultaneously by a microwave interfer- 
ometer and by the integral of the ion loss current, and 
found that ru ~ xo lS m"3. The lifetimes of electrons in- 
jected from the gun were much shorter than the life- 
times of plasma electrons produced by ionization. The 
confinement time of the bulk plasma was about a factor 
of five lower than classical, assuming that the electrons 
were produced by ionization. The injected fast electrons 
left quickly (-0.01 ms), perhaps due to a two-stream in- 
stability, and imparted only about 1-2% of their energy 
to the plasma. The electron current had 0.1-1 MHz 
modulations. Probe measurements showed that the elec- 
tron density profile was not consistent with a theoretical 
profile based on classical transport. '45 

Factors influencing the electron confinement times in- 
cluded the magnetic field geometry, impurities in the 
plasma, ionization of neutral gas and electron mobility 
in the electric field. The effect of the centrifugal force 
on radial electron transport would be less significant in 
large devices. Since the electron mobility is related to 
the neutral gas density, its contribution was significant 
in small devices, but would be reduced in high- 
temperature plasmas, where most of the neutral gas be- 
comes ionized. 1 * 6 

From ion energy distributions measured outside the 
point cusp and ring cusps of the Jupiter i M device at 0 a 
= 6 kV, Azovskii and co-workers : 47 deduced the amount 
of potential sagging in the anode gaps, and found A0 
(ring cusp) = 0.4 kV, A0 (point cusp) = 0.6 kV. These 
values correspond to uncompensated anode electron 
densities n^ a - io '7 m" 3 , while the central electron density 
n. „ - 10 18 m"3. The electronic space charge in the anode 
gaps was partially compensated by ions flowing out the 
cusps. Although the ion potential barrier 0, was lower in 
the point cusp, the ion loss rate there was lower than in 
the ring cusp. From gridded analyzer data they found 
that the ion energy distribution was nearly Maxwellian 
with T, = 40-60 eV T, was slightly higher at the plasma 
edge, where the electric field was large. 1 * 8 

Azovskii and colleagues^ measured the electrostatic 
potential and electron temperature in JupitenM. At the 
end of electron beam injection the potential sag. in the 
point cusp cathode A0 - 1 kV The electron temperature 


reached 25 eV, then dropped rapidly after electron beam 
injection cut-off. Using an 8 keV Ar probe beam, they 
measured the energy spectra of Ar 2+ ions emerging from 
the ring cusp and deduced the time variation of the cen- 
tral plasma potential. When energetic electrons were 
injected into low-pressure hydrogen, the plasma poten- 
tial ( 0i + A0) initially was very large. As ionization built 
up the plasma density, the plasma potential decreased to 
a plateau value on the order of o.x 0a. After the electron 
beam injection was turned off. the plasma density and 
potential decayed to zero in a few milliseconds. 1 * 0 

3.3 Magnetic Multiples 

Nakamura etafr 1 observed enhancement of ion con- 
finement in a multipole magnetic field by a negative 
electrostatic potential well produced by electron injec- 
tion into the plasma. Hershkowitz etal ** 2 observed simi- 
lar effects. 

Lavrent’ev and colleagues 1 ** 1 ** studied the plasma pro- 
duced by injection of 100 mA of 1 keV electrons into a 
pulsed multipole trap. They measured ru ~ io 1 * m"3, con- 
finement time = 50 ps and T» * 300 eV. Using the at- 
tenuation of a calcium ion beam during a 30 ps period of 
electron injection, they observed, a rapid density 
build-up to 4 x 10 [ 7 m" 3 , followed by a decay after elec- 
tron beam cut-off. With probe measurements they found 
a uniform central plasma density dropping off steeply in 
the boundary magnetic field layer. 1 ** Lavrent’ev etaft 6 
found that the plasma lifetime in a multipole electro- 
magnetic trap depended very strongly on the magnetic 
field in the gaps, decreasing from 1.3 ms to 0.1 ms as the 
field was reduced from 0.3 T to 0.15 T. 

Moir etal 1 *7 proposed a toroidal quadrupole MEPC ex- 
periment using soft iron to enhance the magnetic field 
shape and to minimize the required coil currents. 

Maffei etal : * 8 constructed a spherically convergent ion 
focus (SCIF) device based on Bussard’s Polywell concept 
(a polyhedral set of point cusps), with planned values of 
B = 0.2 T, 0 = 20 kV, I = 75 A from three electron guns 
and plasma production by electron cyclotron heating. 

3.4 ATOLL 

The ATOLL toroidal quadrupole magnetic cusp device 
is shown in Figure 21. Typical parameter values are B a = 
0.8T, 0a = 2k V, p = o. 2-0.5 m P a hydrogen, impurity frac- 
tion = 1-2% (dominant mass = 28), ru = 3 x xo '7 m' 3 , 7 J, = 
xoo eV and average ion energy = 100 eV. Ioffe etafo 9 
checked the alignment of the magnetic field with the 
anodes by observing the flow of electrons from an elec- 
tron gun placed in the plasma region. The alignment 
accuracy was better than 0.2 mm. Diagnostics included 
an 8 mm microwave interferometer, currents drawn by 
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electrodes, a charged particle energy analyzer, Langmuir 
probes (including a tiny probe suspended by a 12 pm 
tungsten fiber), a charge exchange neutral atom energy 
analyzer, a pulsed vacuum gauge and a partial pressure 
analyzer. 



In one mode of operation the plasma was sustained by 
electron injection from a gun plus secondary emission at 
the reflecting cathodes. This mode resembled a Penning 
discharge with the end walls outside the cusps. 

In another mode of operation Ioffe and co-workers' 60 
minimized the plasma impurity content by careful wall 
preparation and by the use of very pure hydrogen. Ioni- 
zation of neutral gas made up for particle losses, and the 
electrons were heated by the intense electric field in the 
boundary layer. At high magnetic fields and low pressure 
the discharge could become self-sustaining without elec- 
tron injection or secondary emission from the cathodes. 
This regime had similar plasma parameters, but less 
noise, better ion confinement and higher 7 V, than the 
first regime. This ‘emission-free’ discharge was extin- 
guished by adding 5% of an impurity, such as nitrogen. 
Impurities increased radiation losses and diffusion rates, 
and reduced the potential well depth. 

The plasma potential increased linearly with magnetic 
field and had a nearly parabolic radial distribution. Ions 
produced by ionization at the edge were heated as they 
fell into the potential well, and the average ion energy 
increased with increasing magnetic field. From the 
plasma confinement time it was inferred that the poten- 
tial barrier 0 , impeding ion loss through the cusps was 
very small. Ioffe and colleagues 161 found that the electro- 
static potential in the cusp gaps was equal to the poten- 
tial at the plasma centre, to within experimental accu- 
racy (10%). The hypothesis of a narrow ‘neutral corridor’ 
through the cusps with no confining potential barrier (<P, 
= o) was also consistent with the measured variation of 
plasma potential with magnetic field. 


The electron potential barrier 10 TV,, so electron 
losses were predominantly by cross-field diffusion. From 
the ionization rate Ioffe and co-workers l6 v 6 3 inferred a 
plasma confinement time of 0.18 ms. This was close to 
the theoretical confinement time, based on the ion 
acoustic instability, of 0.12 ms, and much shorter than 
the classical confinement time of 50-150 ms. Similarly, 
the measured electron density profile was close to that 
predicted above for the ion acoustic instability, with a 
gradient scale length on the order of the ion sound Lar- 
mor radius, and far from the profile based on classical 
diffusion. Piterskii and co-workers' 6 ^ 1 ^ observed oscilla- 
tions at 40-80 MHz propagating in the electron dia- 
magnetic drift direction at about io5 m s' 1 , which ap- 
peared to be consistent with ion acoustic instability the- 
ory. According to Lavrent’ev 166 the cross-field transport 
resulting from the observed profiles and oscillations in 
ATOLL should be inwards. 

Later measurements revealed that three instabilities 
were occurring simultaneously in ATOLL diocotron 
oscillations in the cusp gaps, a long-wavelength ion 
acoustic instability in the plasma interior and a lower 
hybrid drift instability at the plasma edge. 

Piterskii etal l6 1 observed diocotron oscillations in the 
cusp gaps at frequencies of 40-70 MHz, with wave- 
lengths of 1-2 cm and peak amplitudes -100V. The wave 
phase velocity was somewhat lower than the electron E 
xB drift velocity in the anode gaps, by a factor of 1.5-2; 
but this difference could be partly due to simplifying 
assumptions used in the analysis, such as the shapes of 
the electron density distribution and electric field dis- 
tribution in the anode gaps. The theoretical maximum 
growth rate of the diocotron oscillations occurs at X, - 
16a, which is within the range of wavelengths observed. 
The instability became stronger with increasing plasma 
potential and decreasing magnetic field, as expected for 
diocotron oscillations. In the central plasma the ampli- 
tude of the diocotron oscillations decayed with distance 
from the anode-gaps with a scale length of about 5 cm, 
as expected theoretically. 

Piterskii and co-workers' 68 , 169 found that the lower hy- 
brid drift instability is dominant in the boundary region 
where the density is low The theoretical values of most 
unstable wavelength (1 cm), frequency (12 MHz) and 
phase velocity (1.2 x ro' m s 1 ) are consistent with ex- 
perimental data. This instability causes anomalous elec- 
tron transport, flattening the electron density profile 
until the density gradient scale length (hl,/VhA exceeds 
the local ion Larmor radius pi. Thus, the confinement 
times in ATOLL were much shorter than classical con- 
finement times, and the anomalous transport was well 
understood in terms of the lower hybrid drift and long- 
wavelength ion acoustic instabilities. 
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3.5 Jupiter-2M 

The proposed Jupiter 2 experiment > 7 ° was to be a linear 
set of ring cusps with B^= 5T (point cusps), ®\ = 25 kV, V 
= 0.5 m 3 , tl, > xo 1 ? m 3, Ti - 1 keV and Te~ ioo ms. Due to 
financial limitations, a one-third scale model of Jupiter 2, 
called Jupiter 2Mm has been built first (Figure 22), with a 
primary goal of studying electron transport. 



Figure 22 — The Jupiter-2M device (left) and electrode region detail (right). 171 

The Jupiter 2M device has seven ring cusps with R = 
0.215 m i n the central region flanked by one small ring 
cusp and point cusp at each end. The axial length be- 
tween the point cusps is 1.30 m. The coils are pulsed 
from a capacitor bank and reach peak fields in about 33 
ms. The .maximum attainable ring cusp and point cusp 
magnetic fields are about 1.94 and 3.88 T, respectively, 
but lower fields have been used in experiments. The 
shape of the vacuum magnetic field is illustrated in Fig- 
ure 23. n 



The plugging electrode positions were adjusted by 
screws until comb probes at three azimuthal locations 
120° apart showed no asymmetry of the electron loss 
currents emerging from the cusps. The total plasma vol- 
ume is 501, and the volume of the uniform central region 
is 161. 


The classical electron diffusional loss rate is 1 ? 2 . m 
J diff = j^D, |\ + 0.5 y~ j + D ei j.F« 0 /? 2 electrons / 5 

( 55 ) 

where D a is the electron diffusion coefficient due to col- 
lisions with neutral atoms, ® p = ®i + A®, D e , is the elec- 
tron diffusion coefficient due to collisions with ions, ru„ 
is the central electron density, T eo is the central electron 
temperature, and F is the geometrical factor having di- 
mensions of reciprocal length. The ® p term represents 
electron mobility in the strong edge electric field. The 
electron confinement time is 


I dW I diff (56) 


where Nt, is the total number of electrons in the plasma. 
For the geometry of the Jupiter2M experiment F = 900/ 


Preliminary low-field operation at B - 0.5 T (ring cusp) 
with a 1.3 kV, 0.12 A electron beam injection into 2 x io' 6 
torr (FI 2 ) produced a plasma with rc = 5 x io 1 / m _ 3, TV, = 13 
eV, Te, > 2 ms. With ®a - 2 kV, the plasma potential was 
about -300 V relative to the grounded chamber walls. 
About 80% of the injected electron beam energy ap- 
peared in the plasma. Although the point cusp anode 
diameters were 25 mm, the ion loss through the point 
cusps was only 6% of the total ion loss rate. 





High-frequency oscillations were not observed in 
Jupiter2M. The large uniform central plasma may serve 
as a buffer to damp out surface oscillations. Vdovin eL 


Plasma Physics and Controlled Fusion 36 (1994), pp. 1539*1593, doi:io. 1088/0741-3335/36/10/001 


28 


al , 1 74 found that the ratio of experimental electron loss 
current to theoretical loss current decreases with in- 
creasing plasma density, as shown in Figure 24. At higher 
plasma densities, the experimental loss rate was about 
twice the classical value. The measured ion loss rate cor- 
responds to an ion potential barrier ®,/T, - 4. 

Auxiliary electrodes were used to suppress secondary 
electron emission from the plugging cathodes. Dia- 
phragms limiting the plasma boundary in the point 
cusps and between the ring cusps served two functions: 

(i) to reduce particle flow onto the anodes, where 
gas emission stimulated by particle bombard- 
ment would tend to fill the gap with cold 
plasma. 

(ii) to study the influence of the diaphragm size on 
plasma transport. 

Ring diaphragms with diameters of 21, 24 and 28 cm 
corresponded to ring cusp diffusion layer widths a = 0.9, 

1.6 and 2.0 mm. 

The theoretical electron density profile depends little 
on the temperature profile. The electron density profile 
is flat in the centre and falls nearly linearly to zero near 
the plasma edge. There is a uniform central plasma with 
radius -5 cm and length -80 cm. As the diffusion layer 
width a is increased from 0.9 to 2.0 mm, the electron 
confinement improves, but ion confinement degrades, 
due to the increase of A (P. There is also an increased 
probability of arcing with large gap widths, due to anode 
bombardment and outgassing. At the optimum width 
a 0 pt, = x .6 mm the maximum plasma density n, = 10' 8 mf 3 
is attained. The electron confinement time deduced 
from the total number of electrons in the plasma di- 
vided by the total electron diffusion loss current is V = 3 
ms. Future plans are to increase the magnetic field to 1-2 
T, and to inject hot plasma with a plasma gun. 

3.6 Impurities 

Impurity effects in MEPC devices have been studied 
very little. Lavrent’ev^ showed the theoretical effect of 
increasing Z e g on the energy confinement time. How- 
ever, nobody has computed impurity source terms, 
transport and build-up. 

A continuous impurity accumulation was observed in 
the ATOLL device during the injection pulse. A typical 
impurity level was 1-2% with the dominant species at 
mass 28. '59 The addition of a few per cent impurity ex- 
tinguished the ‘emission free’ regime of ATOLL. 160 
There have been almost no studies of impurity profiles, 
impurity transport, radiative power loss or methods of 
impurity removal from an MEPC plasma. Impurity stud- 
ies need greater emphasis in the future. As Gormezano39 
stated, “any electrostatic system will constitute a very 


good containment system for heavy impurities. The 
presence of electrodes in environments with high elec- 
tric fields and high particle fluxes may constitute a very 
large source of impurities. This might ultimately be the 
greatest problem.” 

3.7 Comparison of Experiments 

Key parameters of various MEPC experiments are 
summarized in Table 2. In some cases the effective tem- 
perature is taken to be two-thirds of the average particle 
energy The early experiments at Kharkov (S-r, S-3, S-3M 
and S-4M), Quebec, Missouri and Jutphaas were small 
spindle cusp plasmas with weak fields, and neutral gas 
effects were dominant. Such experiments are of less in- 
terest now, but they served to demonstrate the princi- 
ples of MEPC. The multiple point cusp SCIF device was 
shut down by budget cuts before reaching its projected 
operating parameters, but it may be restarted at another 
laboratory 

The JupiterrA device achieved good ion and electron 
temperatures (0.4-1.5 keV). It also demonstrated the 
effectiveness of electron cyclotron heating, although at 
low powers. The JupiteriM experiment was not able to 
fully utilize its high magnetic fields, because the dura- 
tion of the peak field was too short to permit accumula- 
tion of high plasma density Due to their low ratios of 
plasma volume to surface area VI SR < 0.03, spindle cusps 
are not worthwhile considering for future experiments. 

The linear multipole devices demonstrated confinement 
in that geometry and validated theoretical predictions 
of critical loss angles. Due to end losses, these devices 
did not achieve high plasma pressures. 

The ATOLL totoridal quadrupole was carefully built 
and operated, but the electron energy transport was 
anomalous. Detailed studies related the anomalous 
transport to the lower hybrid drift instability (plasma 
edge) and to the ion acoustic instability (plasma inte- 
rior). ATOLL was shut down in the 1980s. 

The Jupiter 2M experiment has attained the most prom- 
ising results. Its demonstration of nearclassical trans- 
port is stimulating interest in MEPC research. 

3.7.1 Density, Temperature 
and Confinement Time Scalings 

The plasma density in spindle cusps was typically io ! 7 to 
10 18 nr 3. The attainment of higher plasma density has 
been a major concern for the reactor prospects of 
MEPC devices. The plasma density appears to scale 
proportional to B 2 , 122 provided that the magnetic field 
duration is long enough for density accumulation. The 
central plasma density it, is found to be significantly 
higher than the electron density n, a in the anode region, 
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but more definitive measurements of the ratio are 
needed at high plasma temperatures. 

The measured electron energy distributions resemble 
shifted Maxwellians, as expected. Computer models in- 
dicate that, if neutral gas effects are not dominant, the 
plasma temperatures scale proportional to the applied 
voltage. Measured ion temperatures have scaled linearly 
with 0a. For the plasma density and temperature in 
MEPC devices, we assume the following scaling 

n ~ Bl ( 57 ) 

T.~T,~ 005 t A(sg) 

The scalings were used to estimate the hypothetical re- 
actor parameters in Tables i and 2. That the tempera- 
tures in Jupiter 2M were below this assumed scaling may 
be due to the short duration of the peak magnetic field, 
the weak heating power and neutral gas effects. 


Table 2A — Parameters of some MEPC Experiments, (italicized values are estimates) 



For small spindle cusp devices the confinement times 
have been t - 0.1 to 1 ms. Little distinction has been 
made experimentally between particle confinement and 
energy confinement. For ring cusps the confinement 
time is assumed to scale roughly in accordance with 
equation (40) 

z diIf ^0.09^~ 

PbPo (40) 

This equation predicts Xdtff= 0.8 ms for Jupiter2M pa- 
rameters of Table 2. The experimentally measured con- 
finement time was 2-3 ms. (If the central electron tem- 
perature were higher than the 15 eV reported from a 
Langmuir probe under the magnet coils, then the theo- 
retical would be longer.) A more precise computation 
specific to the Jupiten2M geometry^ predicts classical 
confinement times about twice the measured value. 
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Table 2B — Parameters of some MEPC Experiments, (italicized values are estimates) 



Reference 

R a 

Be 

Oa 

n (m-3) 

Te 

Ti 

T 

alpe 

P* 

Ve T/R 



(cm) (mm) 

(T) 

(kV) 


(keV) 

(keV) 





Linear multipole, m = 4, Kharkov 

75 

5 

t -5 

0.15 

t.o 








Linear multipole, m = 8, Kharkov 

r 53 > i 54 

6 

1 

t. 

r.o 

*** 

0,0 

0.30 

5.oxxo- 5 

x6 

0.002 

X.XXX 04 

Toroidal quadrupole, ATOLL 

158. 154 

20 

1 

t -5 

,7 

xro. 

0.03 

o.o 7 

X.8XX0. 

81 

O.OOX 

5.3 X 103 

Multiple point cusps, SCIF 

104, 158 



0, 

to 








Linear set of ring cusps, Jupiter 
2M, Kharkov 

173 

22 

r.6 

0.6 

2-5 

xxo. 8 

o.ot 5 

0.03 

3 .0XX0X 

74 

0.005 

3.1 X 104 

Hypothetical reactor, Linear set 
of ring cusps 


4 oo 

2 

6.0 

400 

XI ° 2 ° 

20 

20 

6 

*5 

4-5 

X. 7 XXOS 


3.7.2 Dimensionless Parameters 


Three dimensionless parameter ratios are listed in Table 
2: alp and r v e ,xlR. The a/p^ represents the number of 
electron Larmor radii within one-half of the ring cusp 
anode gap. In MEPC experiments this number has 
ranged from 5 to 100, and for a reactor it would proba- 
bly be in the order of 20-30. The optimum anode gap 
width is of the order of 1.5 to 2 mm, and the electron 
Larmor radius is fixed by the magnetic field and electron 
temperature in a reactor to values of the order of 30-100 
pm. 

The second dimensionless parameter is the ratio of av- 
erage plasma pressure to the peak magnetic field pres- 
sure at the coil where | 3 C is the peak magnetic field at the 
coil = P a - 


B c = 2fi 0 nk 


(T' + T t ) 
B] 


(59) 


(In magnetic confinement fusion it is more common to 
define beta relative to the vacuum magnetic field at the 
plasma centre, but that central magnetic field is zero in 
magnetic cusps, which would yield an infinite value of 
beta.) The values of p c in past MEPC experiments have 
been very low, of the order of 0.001 to 0.1 %. Using the 
density and temperature scalings of equations (57) and 
(58) results in fk = 4.5% for the hypothetical reactor of 
table 2. For comparison, the estimated value for the 
ITER tokamak (Troyon coefficient g = 3 , 1 = 24 mA, mi- 
nor radius = 3 m, B = 6T, B c = 13 T) is |\ = 0.85%. 


The third dimensionless parameter in table 2 is r v e ,x IR, 
which is the ratio of the confinement time to the ap- 
proximate electron transit time through the plasma. 


Present experiments attain confinement values of about 
io4 transit times. Without effective plugging, the plasma 
would only be confined for tens of transit times. 

This number of transit times is relevant to a mode of 
free plasma flow out cusp gaps. It is less relevant to 
MEPC devices, because ions are confined electrostati- 
cally, and their confinement times are roughly equal to 
the electron cross-field diffusion time. Nevertheless, it is 
instructive to show that if the ions were freely flowing 
out cusp gaps, those gaps would be very narrow. Accord- 
ing to the free-flow model, the ion confinement time is 


T = 


r% 

uR8 


( 6 ) 


For Jupiter2M, R = 0.22m, r p = 0.1m, L, = 0.105 m , u ~ 7.6 
x io4 ms -1 and x = 3 ms. Using this equation we estimate 
the effective cusp gap half-width 5 = 2.1 x iot m. Since 
the electron Larmor radius p a = 2.2 x io "5 m, it is seen 
that 6 = p a . If free ion flow were occurring, if would be 
through a very narrow gap. 


3.7.3 Anomalous Transport 

The diocotron instability in the ring cusp regions is be- 
lieved to play a beneficial role by removing anode- 
trapped cold electrons, which tend to screen out the 
anode potential, but other instabilities have deleterious 
effects. The ATOLL toroidal quadrupole cusp experi- 
ment had anomalous transport (hundreds of times faster 
than classical) in all operating regimes, associated with 
the long-wavelength ion acoustic instability and the 
lower hybrid drift instability. The Jupiter 2M multiple 
ring cusp experiment, however, has demonstrated nearly 
classical electron transport. So far, the reason for the 
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difference is not understood. There are several hypothe- 
ses, as follows: 

(i) In Jupiter 2M directed electron injection from 
thermoelectric cathodes was used. According to 
Lavrent’ev, 1 electron injection in ATOLL by 
secondary emission from the cathodes intro- 
duced electrons into the diffusional layer (rather 
than into the central Brillouin region), which 
could foster the two-stream instability. The rate 
of electron production was two orders of mag- 
nitude greater than the electron loss rate by 
classical diffusion, so classical diffusion could 
not accommodate the required outflow rate of 
spent electrons from the trap. 

(ii) ATOLL has two-dimensional curvature and 
only toroidal symmetry in the anode regions, 
while JupitemM has one-dimensional curvature 
and (since h » a) approximately two degrees of 
symmetry. This lack of symmetry in ATOLL 
might be more conducive to growth of the 
harmful instabilities, such as the lower hybrid 
drift waves; or it might somehow tend to reduce 
the beneficial effects of diocotron oscillations 
on cold electrons within the anodes. 

(iii) The large volume of field-free plasma in Jupiter 
2M may serve to damp the growth of instabili- 

More study will be required to falsify or confirm these 
hypotheses. 

3.7.4 Assessment 

MEPC experiments typically have had small radii R - 10 
to 20 cm and plasma volume < 11, except for Jupiter2M 
(501) and the SCIF device. The anode gap half-width 
ranges from 1 to 2 mm. A larger gap could permit exces- 
sive voltage sag A 0 , and a smaller gap would decrease 
the electron diffusion loss time. The magnetic fields 
have generally been below 1 T, except for a few pulsed 
experiments. The applied voltage has been below 10 kY 

Many excellent measurements have been performed in 
MEPC devices. Detailed probe measurements have 
shown that the electron density profile in the ring cusp 
anode has a half-width about one Larmor radius. The 
electrostatic potential profile has been mapped out with 
pulsed ion beam transmission measurements and grid- 
ded energy analyzers. A deep negative potential well has 
been sustained for many milliseconds. The electron ve- 
locity distribution resembles a shifted Maxwellian, as 
expected. These measurements verify the basic theory 
of MEPC. 

The scalings of density, temperature and confinement 
time have been established within the current ranges of 


parameters. The confinement in some experiments is 
nearly classical, when neutral gas effects are taken into 
account. The diocotron oscillations in the ring cusps 
have been thoroughly studied and are well understood. 
The plasma waves associated with anomalous transport 
in ATOLL have been correlated with theoretical predic- 
tions. 

Why some devices have anomalous transport and others 
do not is not yet understood. More measurements are 
needed of the density and temperature profiles, impurity 
distribution and transport, neutral gas distribution, ra- 
diation losses and heat propagation. Two-dimensional 
computer simulations are needed to provide a link be- 
tween analytic theory and experimental data. 

A common thread of deficiencies links almost all the 
past MEPC experiments: either weak or short-duration 
magnetic fields; relatively low voltages (<10 kV); weak 
plasma heating; incomplete diagnostics; inadequate 
computer modeling. These deficiencies have been 
caused by low budgets, not by failures, of the experi- 
menters. 

Table 3 compares the parameters of Jupiter 2M with the 
proposed Jupiter 2 and JupiterT devices. If successful, 
JupiterT could demonstrate the attainment of thermo- 
nuclear parameters in a modesrsized device. 
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Various cusp geometries suitable for electrostatic plug- 
ging were illustrated in Figure 8. Mirror geometries, 
such as a yin-yang field, are also possible for electrostatic 
plugging, but the plasma volume inside the anode- 
bounded flux surface is smaller than for cusp systems. 

4.1 Reactor Design Studies 

Table 3 — Comparison of Jupiter-2M, Jupiter-2, and Jupiter-T experiments. 174 


Parameter 

Jupiter-2M 

Jupiter-2 

Jupiter-T 

Ring cups radius, R (m) 

O.215 

0.67 

1.4 

Fength between ends, L (m) 

1.4 

3-9 

6.4 

Ring cusp magnetic field, B a (T) 

0.5 

2-5 

4.0 

Point cusp magnetic field, B^(T) 

X 

5 

8 

Applied voltage, 0\ (kV) 

2-5 

25 

100 

Electron injection current, A (A) 

0.2 

5 

30 

Plasma density, n, (m‘3) 

0.8 x 10 18 

3 X IO 1 ? 

1020 

Ion temperature, T (keV) 

0.03 

X 

IO 

Plasma volume, V (m 3 ) 

0.05 

0.5 

4 

Electron diffusion current, I jiff (A) 

1.6 

14 

120 

Electron lifetime, e (s) 

0.002 

0.36 

x -3 


4. Reactor 
Concepts 


Laurent’ev'75 considered several reactor geometries for 
MEPC reactors: a sphere, a cylinder and a torus. They 
would typically operate atB a = j T, 0 a = 200 kV, n, - 10 20 
m _ 3, Ti - 30 keV, V - 800 m 3 and power density -4 MW 
mi. They would provide intense sources of neutrons, 
which could be used for production of tritium or pluto- 
nium. 

Pastukhov and Yushmanov'? 6 .^ estimated the power 
gain of an MEPC reactor. At each value of applied po- 
tential 0 a there is a plasma temperature that maximizes 
Q. The optimum gap width makes A 0l0 a = 0.25. When 
the gap is wider, A0 becomes large, reducing the plug- 
ging potentials and 0,. When the gap is narrower, 
electron diffusion losses are too great. At the optimum 
gap width, the ratio of fusion power to input power for ; 
ring cusp plasma reduces to 



where 0 a is in kV and ru L . is the trapped electron density 
in the anode region- Assuming rL,JrL, a = 0.3 and V/SR = 
0.3 in the example case of table 1, it is found that Q = 10. 
This analysis did not include energy losses due to neutral 
gas and impurities, or the beneficial effect of the alpha- 
particle emission current helping to reduce the required 
cathode current. Alternatively, using the parameters of 
Table 1 with Z e g= 3 in equations (4D1) and (3F14) of ref- 
erence 44, we find Q = 13. 

A spherical cusp, shown in Figure 8, would be advanta- 
geous from a physics standpoints 8 ^ The sphere has 
the best surface-to-volume ratio, and the spherical elec- 
trode arrangement might provide some spherical focus- 
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ing of charged particles. If effective, such focusing could 
increase the central plasma density and decrease the 
required magnetic field in the cusps. However, particle 
scattering by waves and by Coulomb collisions and re- 
flections off the scalloped boundary would tend to di- 
minish the degree of focusing. It would be difficult to 
align the electrodes with the magnetic field fines in a 
spherical reactor. 

A torodial hexapole reactor design had B a = 8T, 0 a = 
30okV, n, = io 2 m"3, T = 17 keV and a neutron wall load- 
ing = 1.6 MW nr 2 .*4 The direct capital costs estimated 
for this reactor were comparable to those for tokamak 
and tandem mirror reactors. 1 * 0 

From the results with Jupiter2M, it appears that ion 
losses from the point cusps are not a: serious problem, 
so a linear set of ring cusps would be feasible for a reac- 
tor. Such a reactor would be superior to either the 
spherical cusp or toroidal hexapole cusp, because of its 
simplicity, lower structural mass and easier maintenance. 
A linear MEPC reactor, as illustrated in Figure 25, could 
have the parameters of table 1, with Q = 10, DT fusion 
power density = 3 MW m~3, neutron wall loading = 
3.6MW m 1 , L = 50m, and total fusion power = 4.2GW 



Radiation collimators and voltage dividers would be 
used to prevent high voltage breakdown, as illustrated in 
Figure 26. 



4.2 Applications 

Besides power production, there are other potential ap- 
plications for MEPC devices. Their good confinement 
of high-Z ions would make them well suited for strip- 
ping and confining heavy ions for accelerators, as sugp 
gested by Stix. ,8l > l82 

Nadler etal 18 3 have suggested that small inertial- 
electrostatic confinement (IEC) devices operating in the 
glow discharge mode at 30-50 kV could become useful 
low-cost neutron sources producing >io 6 DD neutrons/s 
for a variety of applications, including detector calibra- 
tion, diagnostics and training. By increasing the voltage, 
current and size, and using DT fuel, yields of 10 10 
neutrons/s might be attainable. Such neutron sources 
could be stronger than 2 5 2 Cf sources (2.6 year half-life), 
and they would be simpler and cheaper than accelerator- 
based neutron sources. 

Miley etal 18 4 found that IEC device operation in the 
‘star’ discharge mode maximizes the effective grid 
transmission factor for ions. An IEC device using D3He 
fuel could generate 14.7 MeV protons for production of 
positron-emitting radioisotopes that are useful for posi- 
tron emission tomography facilities. Miley etal l8 $ pro- 
posed the use of high-voltage collectors in a ‘Venetian 
blind’ configuration to convert the D-3He reaction 
product kinetic energy directly into electricity 
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Miley et al ,ti6 studied hydrogen production by radiolysis 
of water with 14.7 MeV protons from an IEC device. 

The 3.6 MeV alpha energy would be converted directly 
into electricity. The protons would pass through a thin 
chamber wall and interact with steam flowing through 
the cell, splitting about 13 water molecules per 100 eV of 
energy deposited by the protons. The steam would pass 
through a turbine to generate electricity, and then the 
hydrogen would be removed from the condensate. The 
overall system efficiency could be about 50%. 

Bussard' 8 " suggested the use of an electron-beam-heated 
propellant thruster powered by an electrostatic plasma 
device for space propulsion. Miley etal m suggested the 
use of IEC engines with Venetian blind energy conver- 
tors to power five 68~N arcjet-magneto-plasma-dynamic 
thrusters with a specific impulse of 3000 s. 

In pulsed operation, an IEC reactor with R = 5.25 m and 
mass = 5 t could deliver 12 MW of power to the 
thruster. l8 9 Bussard'9 0 proposed the use of a Polywell 
IEC reactor for space propulsion. 

An MEPC plasma with T = 10' 20 keV would be a good 
target for neutral beam injection at W, = 150 keV with a 
beam ion fusion probability - lo' 2 . 1 ? 1 Such a device could 
serve as a powerful neutron source for isotope produc- 
tion, medical therapy or materials testing. 
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5. Issues of MEPC 

Although some recent confinement results appear prom- 
ising, there are still many issues to be resolved before 
MEPC devices can become fusion reactors, as discussed 
below. 

5.1 Self-Shielding of 
Applied Potentials 

The potential drop in point cusp anodes is very large, 
but point cusp ion losses can be kept low by the cen- 
trifugal and magnetic field gradient forces, as demon- 
strated in JupitenM. At B < 6 T, values of n, a s 1.5 x io T ? 
m'3 are required to keep A0 < iookV in the ring cusps. 
The potential sagging in the anode means that high ap- 
plied voltages are needed to ensure substantial and 0,. 

5.2 Neutral Gas Effects 

Neural gas increases cross field transport rates and pro- 
motes high-voltage breakdown. The plasma tends to act 
as a vacuum pump, and in future large MEPC devices 
the neutral gas effects should be less dominant than in 
small spindle cusps. Nevertheless, there will always be 
substantial recycling of hydrogen isotopes and helium 
off the walls, so shields should be used to prevent exces- 
sive flow of neutral gas to the anode regions, where 
build-up of cold trapped electrons would increase A0. 
High-speed vacuum pumping will be required in the 
electrode region. 

5.3 Electrodes 

Optimized plasma confinement by careful alignment of 
each anode was demonstrated in Jupiter 2M. Radiation 
shields and voltage dividers will be needed to sustain 
voltages 0 a > 20okV To avoid warping, the anodes 
should be actively cooled and made from materials with 
low thermal expansion. Although high-voltage holding 
in radiation fields must still be demonstrated in MEPC 
devices, the PSP-2 rotating plasma experiment 1 ? 2 has 
demonstrated the sustainment of 400 kV applied volt- 
age with hot plasma present. 

5.4 Transport 

The reason why Jupiter2M had good confinement and 
ATOLL had anomalous transport is not yet understood. 
Transport in MEPC devices has been studied much less 
than transport in tokamaks and tandem mirrors, so this 
lack of understanding is not surprising. Several hypothe- 
ses have been suggested: the effect of the electron injec- 


tion method, the effects of curvature and symmetry, and 
the effect of a large volume of field-free plasma. The 
extension of the good confinement to larger devices re- 
mains to be tested experimentally A more powerful ex- 
periment is being designed to test the scaling relations 
at higher values of plasma density, temperature and con- 
finement time. 

5.5 Power Gain Ratio 

The estimated power gain ratio for an MEPC reactor 
with nearly classical fuel ion confinement (but no alpha- 
particle confinement), 0.\ = 400 kV, B a = 6 T, R = 4 m, is 
Q = 10. The power gain ratio could be increased by the 
following means: increasing the applied voltage; increas- 
ing the magnetic field; increasing the plasma volume; 
improving alpha energy confinement. 

One means of improving alpha confinement would be to 
add a modest toroidal field to a toroidal multipole device 
(similar to the Tormac configuration). 1 ?? The toroidal 
field would help to confine alpha-particle energy, facili- 
tating ignition, and the resultant magnetic shear might 
help to stabilize the lower hybrid drift instability in the 
plasma edge layer. 1 ?'? Another possibility might be the 
addition of a weak longitudinal magnetic field under- 
neath the surface cusp field of a linear set of ring cusps. 

5.6 Plasma Impurities 

Since alpha particles are poorly confined, helium ash 
accumulation could be slower than in tokamaks. How- 
ever, sputtered wall atoms which penetrate into the 
plasma core would be well confined. Impurity ion accu- 
mulation would dilute the fuel density, cool the electrons 
by radiation losses and diminish the potential well 
depth. One means for impurity removal might be to use 
the centrifugal force of their rotation in a weak longitu- 
dinal field to expand their orbits so that a pumped lim- 
iter could extract them. Another possibility would be to 
expand the point cusp hole at one end of a linear reactor 
and to flow in fresh fuel from the other end, so that the 
weak pressure gradient gradually moved impurities to 
the hole. If effective impurity removal methods could 
not be found, then periodic shutdowns might be re- 
quired for impurity removal, followed by restarts with 
fresh fuel. 

5.7 Materials Issues 

The wall activation, afterheat, accidental release of ra- 
dioactivity and waste disposal issues of MEPC reactors 
would be comparable to those of tokamaks and other 
magnetic confinement devices. The reactor design 
should provide for rapid replacement of damaged walls 
and electrodes. The desirable first wall neutron power 
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flux would probably be about 2-5 MW nr 2 . Higher values 
could cause heat removal and wall lifetime problems, 
and lower values would raise the cost of electricity The 
anodes should be partially shielded from radiation, and 
the electrodes and walls must be actively cooled. Low- 
activation materials, such as purified vanadium chro- 
mium alloys and SiC, should be used if feasible. 

One possible materials combination for an MEPC reac- 
tor is vanadium alloy structure with a circulating Flibe 
(LiF-BeFz) molten salt blanket and coolant. Use of the 
molten salt would avoid the reliability problems of high- 
pressure coolants and the fire hazard of lithium. A small 
amount of lithium could be used to cool first-wall hot 
spots, if needed. The irradiated vanadium alloy could be 
cooled for tens of years, then remelted to remove the 
volatile long-lived isotopes (mostly 4 2 Ar) and recycled, "w 
With careful choice of materials, an MEPC reactor 
would not require active safety systems to prevent 
prompt off-site fatalities in the event of a serious acci- 
dent. 
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6. Conclusions and 
Recommendations 

(0 Purely electrostatic plasma confinement with 

spatial and temporal focusing may be wroth 
re-examination (section 1.1.4). 

(ii) Due to magnetic reflection and electrostatic 
acceleration effects, the central plasma den- 
sity in an MEPC device may be an order of 
magnitude higher than the electron density in 
the ring cusp region. Cold electrons from 
ionization of neutral gas tend to screen out 
the anode potential, but they are probably 
removed from ring cusps by diocotron oscilla- 
tions (sections 2.5 and 2.6). 

(iii) Although anomalous transport was observed 
in the ATOLL toroidal quadrupole device, the 
plasma loss rate in the Jupiter2M Unear ring 
cusp experiment was only about twice the 
classical rate, improving at higher density The 
reason for the good stabiUty of the Jupiter 2 j \1 
plasma is not yet understood, but might be 
related to its method of electron injection, its 
better symmetry, or its large volume of uni- 
form field-free plasma (section 3.7.3). 

(iv) Based on present experimental results, a linear 
set of ring cusps appears to be the most favor- 
able configuration for an MEPC reactor. 

(v) Several decades of small, cheap, pulsed ex- 
periments have demonstrated the basic prin- 
ciples of MEPC, but these experiments have 
been marred by weak or pulsed magnetic 
fields, low voltages, weak plasma heating, in- 
complete diagnostics and inadequate com- 
puter modeling (section 3.7.4). Further pro- 
gress depends upon increasing R, B and ®\ in 
a steady-state device, so that tl., T and Te can 
approach the thermonuclear regime. A new 
experiment is needed to demonstrate the at- 
tainment of interesting parameters (bl, > 3 x 
io T 9 m 3, T > 5 keV) in a steady-state device 
with intermediate size (R = 2 m, L = xo m) and 
reasonable applied fields (B = 5T, ®a = 100 kV). 

(vi) For an MEPC reactor with near-classical 
transport, values of Q = xo are predicted at 
voltages ®.\ = 40okV (section 4.1). Q could be 
enhanced by partial confinement of alpha par- 
ticles. 

(vii) The Polywell and modified Penning trap con- 
cepts are just beginning to be studied. They 


need more research to determine the attain- 
able degree of focusing and to explore related 
concepts that may lead to attractive fusion 
reactors (sections 2.10, 2.11). 

(viii) Radiation shields and voltage dividers should 
be used to prevent high voltage breakdown. 
Voltage holding at 400 kV has been demon- 
strated in a magnetic mirror geometry by the 
PSP-2 rotating plasma experiment (section 
5 - 3 ). 

(ix) MEPC devices may find applications as (sec- 
tion 4.2): electrical power plants; heavy ion 
sources; neutron sources for materials testing, 
medical therapy and isotope production. 

(x) The main issues of concern for MEPC reactor 
development are: electron transport, plasma 
purity and electrode alignment and voltage 
holding. Future research should concentrate 
on resolving these issues. 
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